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Abstract 


This thesis is a collection of published and submitted papers. Each paper is 
the chapter of the thesis and new approach involves proving a pointwise charac- 
terisation of the vectorial infinity Laplacian in the Chapter 2, proving a geometric 
feature of the p — Harmonic and co — Harmonic maps in the Chapter 3, finding an 
explicit oo — Harmonic functions in the Chapter 4, proving two distinct minimality 
principles for a general supremal first order functionals in the Chapter 5. 


In Chapter 2 we introduce the joint paper with N.Katzourakis, which extends 
the result of [56]. Let n, N € N with 2 C R® open. Given H € C?(0 x RN x RY”), 
we consider the functional 


E.o(u,O) := ess supH(-,u,Du), we Wer(O,R%), OE. (1) 
O 


The associated PDE system which plays the role of Euler-Lagrange equations in 
Lis 


Hp(-,u, Du) D(H(-,u, Du)) = 0, 
H(-,u, Du) [Hp(-, u, Du) (Div (Hp(-,u, Du)) —H,(-, u, Du)) 51). 


where [A]+ := Proj R(A) denotes the orthogonal projection onto the orthogonal 
complement of the range R(A) C R% of alinear map A: R" —> R. Herein we es- 
tablish that generalised solutions to (2) can be characterised as local minimisers of 
(1) for appropriate classes of affine variations of the energy. Generalised solutions 
to (2) are understood as D-solutions, a general framework recently introduced by 
N.Katzourakis in [57, 58}. 


In Chapter 3 we present the joint paper with N.Katzourakis and H.Abugirda in 
which we consider PDE system of vanishing normal projection of the Laplacian 
for C? maps u: R” DO —> R®: 


[Du] ‘Au=0 inQ. 


This system has discontinuous coefficients and geometrically expresses the fact 
that the Laplacian is a vector field tangential to the image of the mapping. It 
arises as a constituent component of the p-Laplace system for all p € [2,00]. For 
p = oo, the oo-Laplace system is the archetypal equation describing extrema of 
supremal functionals in vectorial Calculus of Variations in L°. Herein we show 


lv 


that the image of a solution wu is piecewise affine if either the rank of Du is equal 
to one or n = 2 and wu has additively separated form. As a consequence we obtain 
corresponding flatness results for p-Harmonic maps for p € [2, oo]. 


The aim of the Chapter 4 is to derive new explicit solutions to the oo-Laplace 
equation, the fundamental PDE arising in Calculus of Variations in the space L®. 
These solutions obey certain symmetry conditions and are derived in arbitrary 
dimensions, containing as particular sub-cases the already known classes of two- 
dimensional infinity-harmonic functions. 


Chapter 5 is the joint paper with N.Katzourakis. We discuss two distinct mini- 
mality principles for general supremal first order functionals for maps and charac- 
terise them through solvability of associated second order PDE systems. Specifi- 
cally, we consider Aronsson’s standard notion of absolute minimisers and the con- 
cept of oo-minimal maps introduced more recently by N.Katzourakis. We prove 
that C! absolute minimisers characterise a divergence system with parameters 
probability measures and that C? oo-minimal maps characterise Aronsson’s PDE 
system. Since in the scalar case these different variational concepts coincide, it 
follows that the non-divergence Aronsson’s equation has an equivalent divergence 
counterpart. 
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Chapter 1 


Background and motivations 


1.1 Introduction 


Minimization problems have been studied by many mathematician for a differ- 
ent purposes. Most of their efforts were focused in study the relation between 
minimality conditions and partial differential equations (PDEs). One of the way 
to view the minimality as a variational approach, which is the core idea of Calculus 
of Variations. We introduce some fundamental methods of Calculus of Variations 
to solve possibly non-linear PDE, which for a simplicity we have in the following 
form 

L[u] =0, (1.1.1) 


where Llu] is a given differential operator and u is the unknown. The equation 
(1.1.1) can be characterised as the minimiser of appropriate energy functional E[u] 
such that 


The usefulness of this method that now we can proof existence of extremum points 
of the functional energy E]-] and consequently the solution of (1.1.1). One of 
difficulties of described method that in general the minimiser of the functional 
might not be a classical solution of the PDE and the definition of generalised 
solutions is an issue. The generalised solutions that we are using in this thesis 
are not based either on the viscosity solutions which are playing an enormous 
role in nonlinear first and second order PDEs or on the maximum principle or 
on integration-by-parts which helps to “pass derivatives to the test functions”. 
Instead we utilize a recent method that uses the combination of difference quotients 
and Young measures, for precise definition we refer to Section 2.2 and for more 
information about background material of this generalized solutions we refer to 
(ea, ar, A148, OF, FL, 77). 


This thesis is a collection of papers with researcher’s at the University of Read- 
ing, except Chapter 4, which is a single author paper. In all papers we study some 
problems of Calculus of Variations in L° which are in a nutshell a minimisation 
problems with respect to the supremum norm and their corresponding PDEs with 


the method mentioned above. Section 1.2 and 1.3 of this chapter give more details 
on Calculus of Variations in L® and organization of the thesis respectively. 


1.2 Calculus of Variations in L™ 


Calculus of Variations in L® has a long history started in the 1960s by G. 
Aronsson [4-8]. He considered the following variational problem for the supremal 
functional 


E..(u,O) := ess sup |Du(x)|, ue W°(Q,R), O€ ACR” (1.2.1) 


loc 
rEO 


and introduced appropriate L°- notion of minimisers, namely we say the map 
u € W(Q,R) is an absolute minimiser of (1.2.1) if for all O € © and all 


loc 


o € Wy’~(O,R) we have 


E(u,O) < En(u+ ¢,O). (1.2.2) 
Also Aronsson derived the associated PDE 
Au:= Du@Du:D?u=0 inQ, (1.2.3) 


where “®” is a tensor product and “:” is the Euclidean scalar product. Equation 
(1.2.3) is playing the role of L®- analogue of the Euler — Lagrange equation and 
called oo-Laplacian. In particular case when n = 2 and 2 C R? is open bounded 
domain, u(x, ) = |z|3 —|y|2 is the most well-known Aronsson’s solution of (1.2.3) 


which has a C3 regularity and definition of the Hessian on the axes is an issue. 
. . 4 4, Be pe Nice 

However in [9] it was shown that u(x, y) = |x|3 —|y|3 is an absolute minimiser not 

only in Q but in R?. This phenomena was later justified using viscosity solutions. 


The study of vectorial absolute minimisers (i.e. when u € Wy°(Q,RY) and 
N > 2) started much more recently in early 2010s by N.Katzourakis in papers [49— 
59] where he found a new additional term which completes (1.2.3). The associated 


PDE to the functional 


E(u, O) := ess sup|Du(xr)|, ue We"(9,R%), OE NCR", (1.2.4) 


loc 
xrEO 


is the so called co-Laplacian system, 


Att (Du@ Dut [Du}?[DuJ*@ 1) Du = 0, (1.2.5) 


where X € RN*", [X]+ denotes the orthogonal projection on the orthogonal 
complement of the range of linear map X : R" —> RN: 


LX} = Projpyxys- 


Some difficulties of equation (1.2.5) are: the theory of viscosity solutions does not 


work for mappings, coefficients of full system are discontinuous and solutions are 
need to be C? to make a classical sense. However there is a method to reduce 
the regularity of solutions to C! using measures as parameters. For more details 
about this method we refer to Chapter 5. 


For more details about explicit classical solutions of (1.2.3) we refer to the 
introduction of the paper presented in Chapter 4. 


1.3. Organisation of thesis 


The aim of the thesis is to find new classical solutions, derive necessary and 
sufficient conditions and describe a geometric properties of absolute minimisers. 
We have reached our goal by publishing and submitting papers in different fields of 
nonlinear PDEs. Each paper is presented in this thesis as chapter. Every chapter 
below is explained in the outline below. 


Chapter 2 is the joint paper with N.Katzourakis. The estimated contribution 
is 50%. The paper has been accepted at Proceedings of the Royal Society of 
Edinburgh A (Mathematics). We study a pointwise characterisation of the PDE 
system of vectorial calculus of variations in L°. In this chapter we prove that 
generalized solution to PDE 


Hp(-,u, Du) D(H(-, u, Du)) 
H(-,u, Du) [Hp(-, u, Du)]- (Div (Hel, u, Du)) — H,(-, u, Du)) 


0, 
0, 


can be characterized as local minimisers for appropriate classes of affine variations 
to the following energy 


Eoo(u, O) := ess supH(-,u,Du), ue Wer(9,R%), OEQ, 
O 


where n, N € N with 2 C R” open, H € C?(Q x R® x R%"”) is a given, Hp, H,, H, 
denotes the derivatives of H(xz,7, P) with respect to the respective arguments and 
[A]-+ := Proj rays. After an introduction involving a corollary which is a main 
result for classical solutions, we use Young measures and difference quotients to 
define a generalised solutions to fully nonlinear PDE introduced by N.Katzourakis. 
This generalized solutions are called D-solutions. Then we prove two auxiliary 
lemmas that we will use to prove our main result. Finally we prove our main 
result for D-solutions and corollary. 


Chapter 3 presents the joint paper with N.Katzourakis and H. Abugirda [1]. 
The paper has been accepted at the Rocky Mountain Journal of Mathematics. 
We study rigidity and flatness of the image of certain classes of mappings having 
tangential Laplacian. The estimated percentage contribution is 30%. We start 
with a brief introduction bringing an attention to one of results of the paper [52]. 
As generalisation of this theorem we introduce our first result, let Q C R” be an 


open set and n,N > 1 and u € C?(O, RY) be a solution to the nonlinear system 
[Du]+Au = 0 in Q, satisfying that the rank of its gradient matrix is at most one: 


rk(Du) <1 inQ. 


Then, its image u(Q) is contained in a polygonal line in RY, consisting of an at most 
countable union of affine straight line segments (possibly with self-intersections). 
After we show that this theorem is optimal by giving an example that system can 
not have affine image but only piecewise affine. Then we have next theorem as 
the consequence of our first main result which supplement one of the results in 
the paper [53]. Finally we end introduction by our second result which states, let 
Q C R? be an open set and N > 2 let also u: R” DQ —> R® be a classical 
solution to the nonlinear system [Du]+Au = 0 in Q, having the separated form 
u(x, y) = f(x) — f(y), for some curve f € (W3? C?)(R,R%) and some p > 1. 
Then, the image u(Q) of the solution is contained in an at most countable union 
of affine planes in R%. At the end of the chapter we give proofs of our results. 


Chapter 4 is the single author paper published in Journal of Elliptic and 
Parabolic Equations in June 2018. We study explicit oo- harmonic functions in 
high dimensions. We begin our paper with an introduction contained two main 
results for smooth solutions. First result states, let wu: Q C R*? —> R be a C?(Q) 
separated oo-harmonicfunction of the co-Laplace equation in polar coordinates 


1 
—upuz = 0 


2 1 
2 
Up Urr TS — 9 Urleure nae ug tea — pe 


of the form u(r, @) = f(r)g(@). 


(i) Assume |f(r)| = r4 and |g(0)| = e?°, where A and B are any constants, 
then 
A?— A+ B?=0 


or 


(ii) Assume |f(r)| = r4 and |g(0)| = 1g(80) jel CMa then G satisfies the 
following 


— arctan 10) + AP arctan a, if A27-A>0 
A if A=0 
t+c= a . 
— arctan G(t), if A=1 
Git) A-1 Gii)-VA-2?| 
— arctan = 4 sete in ||. if A27-A<0O, 


where c is any constant, provided RHS is well defined. 
or 


(iii) Assume |g(9)| = e2° and |f(r)| = |f (ro) |e! ae , then © satisfies the 


following 


1 6? (t)+B? | B(t)—5 ie pe 
In - arctan if B*-—->0 
2 ©2(t)—(t)+B? 2, /p2_i / 1? 4 
a ©? (t)+B? oe : : if B21 0 
In|t]-+e= 4 3in ®7(1)—B() +B? | " 2eH)—2? A =A 
1 2(t)4+B? 1 B(t)—5—/ ;—-B? 2 1 
3 In €2()-B)+B2| 1 1 In B)—44 rz ‘ if BY — 4 < 0, 


where c is any constant, provided RHS is well defined. 


Finally second result, let n > 2 and wu: Q C R" —> R be a C?(Q) separated 
oo-harmonic function of the co-Laplace equation 


> DiuDjuDju = 0. 


i,j=l 


Then 
|fi(acs)| = |fix?)| eA) for 1 <i Aj <n 


” Bt) de 
Fi(x5)] = [Fi (29)| I 


and 


where F’; satisfies 


F(t F(t 
i+c= ————_——__,, arctan Mt) iM) 


1 
(bay (Lat Cea eso) 


i#j iAj 49 


In the end we provide proofs for the results and numerical experiments. 


Chapter 5 is the joint paper with N.Katzourakis. The estimated contribution is 
50%. The paper has been accepted in the journal Applied Mathematics and Op- 
timization. In this paper we prove vectorial variational principles in L® and their 
characterisation through PDE systems. We start with an introduction involved 
our two main results. 


First result “Variational Structure of Aronsson’s system” says, let u : R" D 
Q —> R™ be a map in C?(0; RY). Then: 


(1) If uw is a rank-one absolute minimiser for (1.1.1) on Q (Definition 5.1.2(i)), then 
it solves 
Hp(-,u, Du) D(H(-,u, Du)) = 0 on Q. (1.3.1) 


The converse statement is true if in addition H does not depend on 7 € R™ and 
Hp(-, Du) has full rank on 2. 


(II) If w has co-minimal area for (1.1.1) on Q (Definition 5.1.2(ii)), then it solves 


H(-, u, Du) [Hp(-,u, Du)}* (Div (Hp(-,u, Du)) — H,(-,u, Du)) = 0 on 0. (132) 


5 


The converse statement is true if in addition for any x € 2, H(z,-,-) is convex on 
R” x RN xn, 


(III) If uw is oo-minimal map for (1.1.1) on Q, then it solves the (reduced) Aronsson 
system 


Aju := Hp(-,u, Du) D(H(-,u, Du)) 
+ H(-,u, Du) [Hp(-, u, Du)}* (Div (Hp(-, u, Du)) — H,(-,u,Du)) =. 


The converse statement is true if in addition H does not depend on 7 € RY, 
Hp(-, Du) has full rank on 2 and for any x € 2 H(z,-) is convex in RN*". 


Second result “Divergence PDE characterisation of Absolute minimiser” says, 
let u : R? D Q —> R® be a map in C1(Q; RY). Fix also O € O and consider the 
following statements: 


(I) w is a vectorial minimiser of E.(-,O) in C}(O;RY)'. 


(II) We have 


max |Hp(-,u, Du) : Dy + H,(-,u, Du)- 2] > 0, 


Argmax{H(-,u,Du) : O} 


for any w € Cj(O;RY)!. 


(III) For any = € C}(O; RY), there exists a non-empty compact set 
Ky, =K C Argmax{H(-,u, Du) : O} (13:3) 


such that, 


(Hp(-,u, Du) Du + H,(-,u, Du) - v)I =i (1.3.4) 


Then, (I) => (I) => (III). If additionally H(z,-,-) is convex on RY x R‘*” 
for any fixed « € Q, then (III) => (I) and all three statements are equivalent. 
Further, any of the statements above are deducible from the statement: 


(IV) For any Radon probability measure? o € P(Q) satisfying 
supp(a) C Argmax{H(-,u,Du) : O}, (123.5) 


we have 
— div(Hp(-,u, Du)o) + H,(-,u, Du)o = 0, (1.3.6) 


in the dual space (Ci(O; RY))*. 


1 We say u € C}(O;R) if u—g € Cj(O;R%), where Co(O;R%) := {py € C'(R";RY) : 
wv =0 on OO}. 

2A Radon measure is a Borel measure that is finite on all compact sets, outer regular on all 
Borel sets and inner regular on all open sets. See [42] for precise definition. 


Finally, all statement are equivalent if K = Argmax{H(-,u,Du) : O} in (III) (this 
happens for instance when the argmax is a singleton set). 


The result above provides an interesting characterisation of Aronsson’s concept 
of Absolute minimisers in terms of divergence PDE systems with measures as 
parameters. The exact distributional meaning of (1.3.6) is 


[ (Hp(-,u, Du) : Dy + H,(-,u, Du) -) do _9 


O 


for all  € C}(O;R%), where the “:” notation in the PDE symbolises the Euclidean 
(Frobenius) inner product in RY*". 


After that, we give a corollary which is a combination of two results in the 
scalar case and for the classical solutions. Then we prove the maximum-minimum 
principle which generalises a corresponding result from [52] and the remark on 
Danskin’s theorem and some of its consequences. Finally, we prove our main 
results using a lemma and the proof of the lemma ends chapter. 


In Chapter 6 we discuss the conclusions and some future work. 


Chapter 2 


A Pointwise Characterisation of 
the PDE System of Vectorial 
Calculus of Variations in L™ 


2.1 Introduction 


Calculus of Variations is the branch of Analysis which deals with the problem of 
finding and studying extrema of nonlinear functionals defined on certain infinite- 
dimensional topological vector spaces, as well as with describing these extrema 
through appropriate necessary and sufficient conditions. Such problems are called 
variational and are ubiquitous in nature, being also of paramount importance 
for other sciences such as Data Assimilation arising in the Earth sciences and 
Meteorology (see [23, 47|). In most applications, the functional one wishes to 
study models some kind of “energy” or “action”. 


Let H € C?(Q x R% x R”) be a given function, where 2 C R” is an open set 
and n, N € N. One of the most standard particular class of functionals of interest 
in Calculus of Variations has the form of 


E(%,0) = [ H(e.u@),Du(a) da 


defined on differentiable maps (i.e. vectorial functions) u: R” D Q —> RY. In the 
above, R“” denotes the space of N x n matrices wherein the gradient matrix 


Dua) = (Dita(x))°— _ ees 


of such maps is valued. We have also used the symbolisations x = (%1,...,%n)!, 
u = (u,...,un)! and D; = 0/0z2;. Latin indices i, j,k, ... will run in {1,...,n} and 
Greek indices a, 3,7, ... will run in {1,..., NV}, even if the range of summation is not 
explicitly mentioned. The simplest variational problem is to search for minamisers 
u of E, sought in a class @ of differentiable maps u, subject to some kind of 


prescribed boundary condition on OQ to avoid trivial minimisers. This means that 
any putative minimiser u € @, if it exists, should satisfy 


E(u, Q) < E(v,Q), for all v € @ with u=v on ON. 


If such a minimiser exists, then the real function t + E(tv + (1 — t)u) has a 
minimum at t = 0 and should satisfy 


Eu +t(v—u)) = 0. 


d 
dt 
By the chain rule, this leads, at least formally, to the next necessary conditions, 
known as the Euler-Lagrange system of Partial Differential Equations (PDE): 


YD i(He,.(,u, Du)) =H,,(-,u,Du), a=1,...,N. 


In the above, the subscripts Hp,,,H,, denote the partial derivatives of H with 
respect to the respective variables P,; and n.. Further, since the integral is additive 
with respect to the domain on which we integrate, it can be easily seen that if wu 
is a minimiser, then 


E(u,O) < E(v,Q), for allv € @ with u =v on 00,7 


where O € Q, namely © is a compact subset of Q. The above weaker condition 
still suffices to derive the Euler-Lagrange system and any putative u satisfying it 
is called an absolute (or local) minimiser. 


The above discussion, although completely formal, nonetheless captures the 
quintessence of Calculus of Variations. However, one needs to use hardcore analytic 
tools to make rigorous the above formal reasoning. In particular, a central problem 
is that the minimisers are sought in a class of at most once differentiable maps, 
which the PDE is of second order and one has to devise a way to make sense 
of the PDE weakly, since second derivatives of u may not exist! Such objects 
are called generalised solutions. Finding a efficient concept of generalised solution 
which allows one to prove that such a generalised object in fact exists and study its 
properties is a highly nontrivial part of the problem. A particular relevant question 
of great interest is to identify conditions on H allowing to characterise variationally 
the PDE system in terms of the functional, namely to provide sufficient as well as 
necessary conditions. 


In this paper we are interested in the variational characterisation of the PDE 
system arising as the analogue of the Euler-Lagrange equations when one considers 
vectorial minimisation problems for supremal functionals of the form 


E(u, O) := ess supH(z,u(z),Du(z)), O EQ, (22.1) 


rEO 


defined on maps u: R” D OQ —> R%. This is in the spirit of the above discussion, 
but for the modern class of functionals as in (2.1.1). The scalar case N = 1 first 


arose in the work of G. Aronsson in the 1960s |[6, 7] who initiated the area of 
Calculus of Variations in the space L°®. The field is fairly well-developed today 
and the relevant bibliography is vast. For a pedagogical introduction to the topic 
accessible to non-experts, we refer to [54]. 


The study of the vectorial case N > 2 started much more recently and the full 
system (2.1.2)-(2.1.4) first appeared in the paper [49] in the early 2010s and it is 
being studied quite systematically ever since (see [50],[53],[52],[51],[55], [57],[56], as 
well as the joint works of N.Katzourakis with Abugirda, Pryer, Croce and Pisante 
[2, 31, 63, 64]). The appropriate class of maps to place and study the functional 
is the Sobolev space W1(Q, RY) of L® maps with L® derivative defined a.e. on 
Q (see e.g. [37]). The direct extension of the concept of absolute minimisers for 
(2.1.1) reads 


Eo(u,O) < En(u+¢,0), OE, 6€Wy™(O,R) 


and was introduced and studied by Aronsson in the context of the scalar case. The 
subscript nought means that @ = 0 on 0O. The associated PDE system arising 
from (2.1.1) as a necessary condition is 


Fxo(-,u,Du,D?u) = 0 in Q, (2.1.2) 


where 


Fat US eR 2’ SR 


is the Borel measurable map given by 
Foo(,9, P,X) := Hp(x,n, P)(Hp(o,9, P): X + Hy(x,n, P)'P + Ho(x,n, P)) 


+ H(a,n, P) [Hp(2,n, P)| (Hpp(o,n, P) :X + Hp,(2,n, P):P 
(2.1.3) 


+ Hp,(x,n, P) :1- H,(2,7, P)). 


In the above, RX n° symbolises the space of symmetric tensors wherein the hessian 
of u is valued: 


qeeey 


ei 


Further, [A] denotes the orthogonal projection onto the orthogonal complement 
of the range R(A) C R% of a linear map A: R" —> R*: 


[A]}~ := Projacays- (2.1.4) 
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In index form, F,, reads 
Foo(Z, 1), DP, X)a — > Hp,,(2, 1), -) (= Hp (2; 1), P)Xgij = 9 Hy, (2, 1), P) Pai 
a Bj B 


+ HCP) - Han PS. [Hp(z,n, P)laa- 
B 


(= Hp, :P53(t;;P) Xpiy + D | Hpying (a, 7, P) Pas 
tJ i 


= > Nealon, PY) _ Hy, (x, 7, P)) 


where a = 1,..., N. Note that, although H is C?, the coefficient [Hp(-, u, Du)]+ is 
discontinuous at points where the rank of Hp(-,u, Du) changes. Further, because 
of the perpendicularity of Hp and [Hp]+ (that is [Hp]+Hp = 0), the system can 
be decoupled into the two independent systems 


Hp(-,u, Du) D(H(-, u, Du)) 
H(-,u, Du) [Hp(-, u, Du)]- (Div (Hel, u, Du)) — H,(-, u, Du)) 


0, 
0. 


When H(z, 7, P) = |P|? (the Euclidean norm on R“” squared), the system (2.1.2)- 
(2.1.4) simplifies to the so-called co-Laplacian: 


Dygth = (Du@ Dut [Du}?[DuJ*@ 1) Daal (2.1.5) 


In this paper we are interested in the characterisation of appropriately defined 
generalised vectorial solutions u : R" D Q —+ R% to (2.1.2)-(2.1.4) in terms of 
the functional (2.1.1). It is well known even from classical scalar considerations 
for N = 1 that the solutions to (2.1.2)-(2.1.4) in general cannot be expected to be 
smooth. In the scalar case, generalised solutions are understood in the viscosity 
sense (see [28, 30, 54]). Since the viscosity theory does not work for (2.1.2)- 
(2.1.4) when N > 2, we will interpret solutions in the so-called D-sense. This is 
a new concept of generalised solutions for fully nonlinear systems of very general 
applicability recently introduced in [57, 58}. 


Deferring temporarily the details of this new theory of D-solutions, we stress the 
next purely vectorial peculiar occurrence: it is not yet known whether Aronsson’s 
variational notion is appropriate when min{n,N} > 2. In the model case of 
(2.1.5) and for C? solutions, the relevant notion of so-called co-Minimal maps 
allowing to characterise variationally solutions to (2.1.5) in term of u ++ ||Du||z-0;.) 
was introduced in [52]. These findings are compatible with the early vectorial 
observations made in [17, 18], wherein the appropriate L™ quasi-convexity notion 
in the vectorial case is essentially different from its scalar counterpart. In the recent 
paper [56] a new characterisation has been discovered that allows to connect D- 
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solutions of (2.1.5) to local minimisers of u +> ||Du||;-0(.) in terms of certain classes 
of local affine variations. This result offered new insights to the difficult problem 
of establishing connections of (2.1.1) to (2.1.2)-(2.1.4). 


In this paper we generalise the results of [56], characterising general D-solutions 
to (2.1.2)-(2.1.4) in terms of local affine variations of (2.1.1). Our main result 
is Theorem 2.3.2 that follows and asserts that D-solutions to (2.1.2)-(2.1.4) in 
C1(Q, RY) can be characterised variationally in terms of (2.1.1). The a priori C? 
regularity assumed for our putative solutions is slightly higher than the generic 
membership in the space W1°°(Q, R), but as a compensation we impose no con- 
verity of any kind for the hamiltonian H for the derivation of the system. 


In special case of classical solutions, our result reduces to the following corollary 
which shows the geometric nature of our characterisation!: 


2.1.1 Corollary [C? solutions of F,, = 0] 


Let 2 C R” be open, u € C?(, RR”) and H € C?(Q x R® x R*”). Then, 


E(u, O) < Eo(u+A, OQ), 


Fx.(-,u,Du,D?u) = 0 inQ => | i 
VOEQ, VAE (AGP UAG™) (u). 


Here All (u), Ao (u) are sets of affine maps given by 


D?A =0, A(x) = 0 and exist € € RYand 
Alli (u) = <A: R" >R” | x € Olu) st. the image of A is parallel 
to the tangent map of €H(-,u, Du) at x 


? 


D?A = 0 and there exists x € O(u) s.t. the 
As™(u) = 4 A: R" > RN | image of A is normal to Hp(-,u, Du) at x 
and A'Hp(-,u, Du) is divergenceless at x 


O(u) := Argmax {H(-, u, Du) : O}. 


This paper is organised as follows. In Section 2.2 that follows we record all the 
basic facts needed regarding the concept of our D-solutions, namely our notion 
of generalised solution required to make rigorous sense of (2.1.2)-(2.1.4). We also 
include a quick introduction to the analytic setup of so-called Young measures, on 
which D-solutions are based. We also give two simple auxiliary results which are 
utilised in the proof of our variational characterisation. Finally, in Section 2.3 we 
state and prove our main result. 


'We caution the reader that the statement of Corollary 2.1.1 sacrifices precision for the sake 
of clarity. The fully precise statement is that given in the main result, Theorem 2.3.2. 
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2.2 Young measures, D-solutions 
results 


2.2.1 Young Measures 


and auxiliary 


Let Q C R” be open and K a compact subset of some Euclidean space RY The 
set of Young measures Y (Q, K) forms a subset of the unit sphere of a certain L° 
space of measure-valued maps and this provides its useful properties, including 
sequential weak* compactness. More precisely, Y (Q, K) is defined as 


Y (Q,K) := {v >: 0 — A(R) | [v(-)|(“) € L°(Q) for any open U C K}, 


where “(K) is the set of Borel probability measures on K. To see how it arises, 
consider the separable space L* (Q, C (K)) of Bochner integrable maps. This space 
contains Carathéodory functions ® : Q x K —> R (namely functions for which 
(-, X) is measurable for all X € K and ®(z,-) is continuous for a.e. x € 2) which 


satisfy 


| Pll z2(a,ccw)) = f I lew < 00. 


We refer e.g. to [35, 43, 77] and to [56, 57] for background material on these spaces. 


The dual space of this space is LS. (Q, M(R)), namely 


(119, c(K))) = L®(2,.M(K)). 


This dual Banach space consists of Radon measure-valued maps 2 3 2 +> v(x) € M(K) 


which are weakly* measurable, in the sense that for any open set U C K, the func- 
tion x + |y(x)|(Y) is in L°(Q). The norm of the space is given by 


lly ll ze, (@,mcK)) = ess sup ||v(x)||, 
rE 


where “|| - ||” denotes the total variation. It thus follows that 


Y (2,1K) = {v € L&(2,M(K)) : v(x) € A(R), 


for a.e. x E ah. 


2.2.2 Remark [Properties of Young Measures] 


We note the following facts about the set Y (Q, K) (proofs can be found e.g. in 


ii): 


i) It is convex and sequentially compact in the weak* topology induced from L°%.. 


ii) The set of measurable maps V : R” > 2 —> K can be 
of it via the embedding V +> oy, dy (x) := bya): 
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identified with a subset 


iii) Let V’, V® : R" D> Q —> K be measurable maps, i € N. Then, up the passage 
to subsequences, the following equivalence holds true as i + oo: V’ —> V™ aze. 
on 2 if and only if dyi + dye in Y(Q, K). 


2.2.3 D-solutions 


We now give some rudimentary facts about generalised solutions which are required 
for the main result in this paper. For simplicity we will restrict the discussion to 
n = 1 for maps u: R D Q —> RN with 2 an interval. The notion of D-solutions is 
based on the probabilistic interpretation of limits of difference quotients by using 
Young measures. Unlike standard PDE approaches which utilise Young measures 
valued in Euclidean spaces (see e.g. [24, 37, 41, 43, 67, 71, 77]), D-solutions are 
based on Young measures valued in the 1-point compactification RY := R% U {oo} 
(which is isometric to the sphere SY). The motivation of the notion in the case 
of C' solutions to 2nd order fully nonlinear systems is the following: suppose 
temporarily u € C?(Q, RY) is a solution to 


Pate) wey =O, we, (2.2.1) 


where F : Q x RX x R® x RY — RX is continuous. Let D!” be the usual 
difference quotient operator, i.e. D'u(x) := +[v(a +h) — v(a)], c €Q,h £0. It 
follows that 


F(e,w(2), w(x), lim Diu'(z)) =0, wen. (2.2.2) 
— 
Since F is continuous, (2.2.1) is equivalent to 


lim F(x, u(x), w(z),D™u'(2)) =0, «EQ. (223) 
The crucial observation is that the limit in (2.2.3) may exist even if that of (2.2.2) 
does not, whilst (2.2.3) makes sense for merely C! maps. In order to represent 
the limit in a convenient fashion, we need to view u” and the difference quotients 
D1"! as probability-valued maps from Q to F (RY My given by the respective Dirac 
masses © ++ dp2y(z) and © ++ Opiny(e)- The exact definition is as follows: 


2.2.3.1 Definition [Diffuse Hessians] 


Let wu: R° DO — R* be in WE°(Q, RY). Let also Dt" denote the difference 
quotient operator, i.e. D1” := (oo pe Di") and Dp"v = z[v(- + he’) — v], hs 
0. The diffuse hessians Du of u are the subsequential weak* limits of the 
difference quotients of the gradient in the set of sphere-valued Young measures 
along infinitesimal sequences (h,)%°, (i.e. lim hg= 0) 


} —3D*u in ¥(Q,RX"), as k + oo. 


piv, Du 
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The above means for any ® € L1(Q, C(K)), we have 


[ fee Xbb dO di [ [ e@xapeacx) dz, ask — oo. 


Note that the set of Young measures is sequentially weakly* compact hence 
every map as above possesses diffuse 2nd derivatives. 


2.2.3.2 Definition [D-solutions to 2nd order systems] 


Let Q C R” be an open set and F : 2 x R® x RX” x RNn? — R” a Borel 
measurable map which is continuous with respect to the last argument. Consider 
the PDE system 


F (-,u, Du, D’u) =) on (, (2.2.4) 


We say that the locally Lipschitz continuous map u : R” D Q —> R% is a D- 
solution of (2.2.4) when for any diffuse hessian D?u of u, we have 


sup |F (x, u(x), Du(x), Xz)| =0, ae ren. (2.2.5) 
X.€ supp, (D2u(z)) 


Here “supp,” symbolises the reduced support of a probability measure excluding 
infinity, namely supp,(V) := supp(V) \ {co} when 0 € P( Ry. 


We note that D-solutions are readily compatible with strong/classical solutions: 
indeed, by Remark 2.2.2iii), if w happens to be twice weakly differentiable then we 
have D?u(x) = dp2u(2) for a.e. x € Q and the notion reduces to 


sup |F (x, u(x), Du(x), X,)| =0, ae reEQ, 


Xe€ supp(9p2.4() ) 


thus recovering strong/classical solutions because supp(dp2u()) = {D7u(2)}. 


2.2.4 Two auxiliary lemmas 


We now identify two simple technical results which are needed for our main result. 


2.2.4.1 Lemma 


Suppose 2 C R” is open, u € C'(Q, R”) and H € C?(R” x R® x R®”). Fix O € 
and an affine map A: R"” —+ R™. We set 


O(u) = {e EO : H(z, u(x), Du(z)) = Ex(u, o)}. 
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a) If we have E,.(u, O) < E..(u+tA, O) for all t > 0, it follows that 


max { Hp(z,u(z), Du(z)) :DA(z) + H,(z,u(z), Du(z)) A(z) > 0. 


z€O 


66,99 


RNn 


In the above and “-” denote the inner products in and RY respectively. 


b) Let « € O and 0 <« < dist(z,0O). The set 


O.(x) = {y €O : H(y,u(y), Du(y)) < H(2,u(2), Du(2)) + 1) B-() 
(where “(-)°” denotes the interior) is open and compactly contained in O, whilst 
Eoo(u, O-(2)) = H(x, u(x), Du(x)), 


whenever O.(x) 4 0. 
Note: The proof does not use affinity of map A. 


2.2.4.2 Proof of Lemma 2.2.4.1 


a) Since E(u, O) < E,(u+tA, O), by Taylor-expanding H, we have 
0O< mae H(-,u +tA, Du + tDA) _ mex H(., u, Du) 
a {H(-,u, Du) + ¢H,(-,u,Du)-A + tHp(-,u,Du):DA 
+ O(?|A?? + P|DAP) } a H(-, u, Du) 


IA 


t max {Hy (+ Du)-A + Hp(-,u, Du) : pa} + O(t?). 
O 


Consequently, by letting t — 0, we discover the desired inequality. 


Item b) is a direct consequence of the definitions. 


Next, we have the following simple consequence of Danskin’s theorem [34]: 


2.2.4.3 Lemma 


Given an open set 2 C R", consider maps u € C1(0,R%) and H € C?(R” x 
RY x RX”) such that H(z,-,-) is jointly convex for any x € , an affine map 
A:R" — R* and O € 1. We define 


r(X) := Enj(uwt+AA, O) — E(u, O), A>O. 
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Let also O(u) be as in Lemma 2.2.4.1. Then, r is convex, r(0) = 0 and also it 
satisfies 


Dro) = max {Hp(-, 1, Du) ‘DA - Hels, Dy): A}, 


where Dr(at) := lim inf Matar ria) is the lower right Dini derivative of r at a. 
rA>0 


2.2.4.4 Proof of Lemma 2.2.4.3 


The result is deducible from Danskin’s theorem (see [34]) but we prove it directly 
since the 1-sided version above is not given explicitly in the paper. By setting 


RO, y) = H(y, u(y) + AA(y), Duly) + \DA(y)) 


we have r(A) = max,<g R(A, y) — max,.@ R(0,y), whilst for any A > 0 the maxi- 
mum max,.G R(A, y) is realised at (at least one) point y € O. Hence 


(r(A) — r(0)) = x [max R(A,y) — max ROO, ¥)| 


yeO yeO 


|r 


= S[ROW) — RO,¥%)| 


A 
x [ROW — R099) + (RAY?) — R(O.2"))] 


and hence 


— 


5(r) = 110) = <(RO0") = ROW"), 


where y° € O is any point such that R(0,y°) = maxg R(0,-). Hence, we have 


Dr(0*) = lim inf ey — r(0)) 


A>0t 


IV 


max {li inf * (RO, y°) — RO, v")) } 


y2cO | A>Ot 


yEO(u) A>0+ 


= max {limin 5 (720.9) = R(0,y)) } 


= max {lim nt = (H( -,u+AA, Du + ADA) — H(-, U, Du)) \ 
O(u) A30t A 


= max {limin + (1H, u,Du) + A Hy (-; 2; Dit) -A + AHp(-,u,Du): DA 


O(u) A>0+ 


+ O(|ADA|? + |AA]?) — H(u,Du)) 
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and the desired inequality has been established. Finally by convexity of H we have 
for any x,y > 0 and any ¢ € [0,1] 


r(tz+(1—d)y) = E(t(ut+ 2A) + (1—t)(ut yA), O) — E(u, O) 
<tE.(u+ 2A, O)+(1-—t)E.(u+ yA, O) — E(u, O) 
< ir(x) + (1 — #)r(y) 


Let us record the next simple inequality which follows from the definitions of 
lower right Dini derivative, in the case that H(z, -,-) is jointly convex for any x € 2. 
This is 

r(A) — r(0) > Dr(0t) A, (2.2.6) 
for all A > 0. 


2.3 The main result of the Chapter 2 


Now we proceed to the main result of the paper, the variational characterisation 
of D-solutions to the PDE system (2.1.2) in terms of appropriate variations of the 
energy functional (2.1.1). We recall that the Borel mapping F,, : 0 x RY x RX” x 
RN” _. RN is given by (2.1.3)-(2.1.4) and Q C R” is a fixed open set. 


2.3.1 Notational simplifications and perpendicularity con- 
siderations. 


We begin by rewriting F,,(-,u,Du,D?u) = 0 in a more malleable fashion (see 
(2.1.3)). We define the maps 


FaeieP Xs) = Hpp(z, 7, P) oak Hp, (x, 1, P) oF es Hp,(z,7, P) : iF (2.3.1) 


Fili(e,,P,X) := Hp(2,n, P): X + H,(2,9,P)'P + H,(2,n, P) (2.3.2) 
and these are abbreviations of 
FAP, Dae = > Mpapy 7,P)Xpg + Means es P) Pas 
Byt,J Bi 
+ S Hea, P) ; 


Fi (a,n, P,X); = S$ Hp,, (2,7, P)X pig + S| Hy (2,9, P) Pai + He,(2,, P). 
Bj B 
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Note that F3(x,n, P,X) € RY, whilst Fho(ax, n, P, X) € R”. By utilising (2.3.1)- 
(2.3.2), we can now express (2.1.3) as 


Fx DZ; 1); E X) = Hp(2, 7); P) Fil, @, 1); a X) Fe Hg, 7); P) : 
. [Hp(2,9, Pt (FE(@0, PX) _ Hy, (2,7, P)). 
Further, recall that in view of (2.1.4), [Hp(x,, P)]-+ is the projection on the or- 


thogonal complement of R(Hp(x,7, P)). Hence, by the orthogonality of [Hp(x,7,P)]+- 
(FE (#,,P,X) — Hy(x,,P)) and Hp(e,, P) F(a, n,P, X), we have 


Foo(2,7, P,X) =0, for some (2,7, P,X) € Qx RY x RX™x RN, 
if and only if 


Hp(z, 1); P) F(z, 1); Hae X) a 0, 
H(x,9, P) [Hp(@.9, P)Y (F(@,0, P,X) -Hn(2,7, P)) 


I 
oS 


Finally, for the sake of clarity we state and prove our characterisation below only 
in the case of C! solutions, but due to its pointwise nature, the result holds true for 
piecewise C! solutions with obvious adaptations which we refrain from providing. 
We will assume that the Hamiltonian H satisfies 


{Hp(2, n, ‘) = 0} C {H(«,n, ‘) = Oo}, (2,9) EN x RN. (2.3.3) 
We will also suppose that the next set has vanishing measure 
{2 EQ: B,,(z) () {h > h(x)} is dense in B,,(z)} = (), (2.34) 


where r, = dist(z,0Q) and h = H(.,u, Du). This assumption is natural, in the 
sense that it is satisfied by all know examples of explicit solutions (see [55, 63-66]). 
It is trivially satisfied if h has no strict local minima in the domain. 


Lets examine three examples for conditions (2.3.3) and (2.3.4). For all examples 
H(z,7,P) = |P|?. Clearly (2.3.3) is satisfied. Remains to show that actually 
condition (2.3.4) holds for our three examples. 

Example 2.3.1. u(x, y) = |z|3 — |y|3 is well-known explicit solution and let Q = 
[—1,1]?. The function h = |Du|? = 16 (|a°|5 a lyl3) has only one point of local 
minimum at origin which means set {h > h(0)}NB,,(0) is the dense in the ball 
B,,(0). Let a point (x,y) be different from origin then it easy to check that set 
{h > h(x, y)} OB,,,(v, y) ts not a dense in the B,,.,(@,y) and as the result we 
have (2.3.4). 

Example 2.3.2. Let Q = (0.1, 1]? and u(x,y) = \/x?2 + y? is the conic solution. 
The function h = |Du|? = 1 for any point of Q. So clearly (2.3.4) is satisfied. 
Example 2.3.3. Using notation e = (cost,sint) we have vectorial solution 
u(z,y) = e* —e¥ on Q = [-1,1]? which is Eikonal, namely |Du|? = |D,u|? + 
|D,u|? = 2. So clearly (2.3.4) is satisfied. 
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Our main result is as follows: 


2.3.2. Theorem [Variational characterisation of the PDE 
system arising in L™] 


Let 2 C R” be open, u € C'(Q,R) and H € C?(Q x R” x R®”) a function 
satisfying (2.3.3) and suppose that (2.3.4) holds. Then: 


(A) We have 
Foo (-,u, Du, D?u) = 0 in Q, 


in the D-sense, if and only if 
Exo(u,O) < Ex(ut+A,O), VOEQ, VA AL™(u) J AG™(u). 


For the sufficiency of the PDE for the variational problem we require that H(z, -, -) 


||,co 


be convex. In the above, the sets AG*(u),A5™(u) consist, for any O € , by 
affine mappings as follows: 


D?A =0, A(x) = 0 for x € O(u), & exist € E R® 
All (u) = 4A: R" > RY | Du € Y(O,RN”) & X, € supp, (D?u(z)) U 
st.: DA = € @ F(x, u(x), Du(x), X.) 
{A ee ae const} 


and 

D?A = 0 & there exist x € O(u), Du 

€ Y (2, RX’) & X, € supp, (D?u(z)) 
ance, 

s.t.: A(x) € R(Hp(z, u(x), Du(e)) ) 

& DA€ L(x, A(x), Xz) 

{A Res RX A= const } 


As (u) =< A: R" > RN 


where O(u) defined in lemma 2.2.4.1 and Y(x,7, X) is an affine space of N x n 
matrices, defined as 


{@ € RX"|Hp(z,u(z),Du(x)): Q = —n- Fz,(x,u(z),Du(x),X) \. 
Lax) = if Hp(zx, u(x), Du(x)) 4 0, 
{0}, if Hp(x, u(x), Du(x)) = 0, 


for any («,n,X) €Qx RN x RN”, 


(B) In view of the mutual perpendicularity of the two components of F,, (see 
(2.3.1)-(2.3.2)), (A) is a consequence of the following particular results: 


Hp(-,u, Du) Fl (.,u, Du, D?u) = 0 in Q, 
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in the D-sense, if and only if 
E(u, O) < Ex(ut+A,O), VOED, VA AL™(u) 
and also 
H(-,u, Du) [Hp(-,u, Du}! (FE, u, Du, D'u) — Hy(-,u,Du)) = 0 in 
in the D-sense, if and only if 


E(u, O) < E,(u+A,O), VOEQ, VAE AG™(u). 


We note that in the special case of C? solutions, Corollary 2.1.1 describes the 
way that classical solutions u : R” D Q —> RN to (2.1.2)-(2.1.4) are characterised. 


2.3.3 Remark [About pointwise properties of C' D-solutions|] 


Let u: R” D Q —> R* be a D-solution to (2.1.2)-(2.1.4) in C1(Q, RY). By 
Definition 2.2.3.2, this means that for any D?u € Y(Q, Re), 


Fa (x, u(x), Du(x), X.) =0, ae. x EQ), V X, € supp, (D?u(z)). 


By Definition 2.2.3.1, every diffuse hessian of a putative solution is defined a.e. on 
Q. as a weakly* measurable probability valued map R" DQ —> P( RN wy {oo}). 
Let 23 2+ O, € RX” be any selection of elements of the zero level sets 


{X« Re Fol tte), Dua), Xx) = 0}. 


By modifying each diffuse hessian on a Lebesgue nullset and choosing the repre- 
sentative which is redefined as D?u(x) = 60, for a negligible set of x’s, we may 
assume that D?u(zx) exists for all x € Q. Further, given that Du(x) exists for all 
x € Q, by perhaps a further re-definition on a Lebesgue nullset, it follows that u 
is D-solution to (2.1.2)-(2.1.4) if and only if for (any such representative of) any 
diffuse hessian 


Fx.(v,u(x),Du(x),X,) = 0, VreQ, VX, € supp,(D*u(z)). 


Note that at points x € Q for which D?u(x) = 6;..} and hence supp, (D?u(x)) = 0, 
the solution criterion is understood as being trivially satisfied. 


2.3.4 Proof of Theorem 2.3.2 


It suffices to establish only (B), since (A) is a consequence of it. Suppose that 
for any O € 2 and any A € Aa (u) we have E.(u,O) < E.(u+A,O). Fix a 
diffuse hessian D?u € Y(Q, Re) a point z € O such that supp,(D?u(z)) 4 0 
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and an X, € supp,(D?u(zx)). In view of (2.3.1), if Hp(z, u(x), Du(x)) = 0, then, 
by our assumption on the level sets of H, we have ie bec u(x), Du(ax ) = 0 as well 
and as a consequence we readily obtain 


H(v, u(x), Du(2)) [Ho(x, (2), Du(e)) 


: (F(z, u(x), Du(e), X.) — Hy (x, u(x), Du(x))) = 0 (2.3.5) 


is clearly satisfied at x. If Hp(x,u(x),Du(x)) 4 0, then we select any direction 
normal to the range of Hp(x, u(x), Du(x)) € RN”, that is 


Ny € R(Hp(2,u(x),Du(x))) C RN 


which means nj Hp(a, u(), Du(z)) = 0. Of course it may happen that the linear 
map Hp(a,u(z),Du(z)) : R" —> RY” is surjective and then only the trivial 
nz = 0 exists. In such an event, the equality (2.3.5) above is satisfied at x because 
[Hp(2, u(x), Du(x))]> = 0. Hence, we may assume n, # 0. Further, fix any matrix 
N, in the affine space #(x,nz, Xz) C R®”. By the definition of Y(x,nz, Xz), we 
have 


Hp( x, tz); Dult)) i Ng = ey Fn (2, ula); Duta) 
Consider the affine map defined by 


A(z) := ng + N,(z-—2), z ER”. 


We remark that tA € Aj~°(u) for any t € R. Indeed, this is a consequence of our 
choices and the next homogeneity property of the space f(x, 17, X): 


(et RX) Ht ei Xx), CER 


Let ¢ > 0 be small, fix x € ( and let us choose as O the domain O,(x) defined in 
Lemma 2.2.4.1b). Our assumption (2.3.4) implies that O.(a) # @ for ae. x € 1. 
In view of the above considerations, we have 


Eco (u, O2(2)) < Eoo(u + tA, O-(2)). 
By applying Lemma 2.2.4.1a), we have 


Os pes {Hp(z,u(z), Du(z)) : DA(z) + H,(z, u(z), Du(z)) - A(z)} 


avi Hp(z, u(x), Du(x)) : Nz + H,(2, u(x), Du(x)) + ne 
= —Nz- t= (x, u(x), Du(x),X,) —H, (x, u(x), Du(e))). 


As a result, we have 


Ny * Ge (x, u(x), Du(x), X,) - H,(2, w(x), Du(z)) J <0 
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for any direction n, L R(Hp (uz), Du()) ) and by the arbitrariness of n,, we 
deduce that 


[Hp(2, u(x), Du(2)) + (FE (w, u(2), Du(2), Xe) — H,,(x, w(e), Du()) ) =i), 


for any D?u € Y(Q, nN), a €Q and X, € supp,(D?u(z)), as desired. 


For the tangential poe of the system we argue reap Suppose that 
for any O € Q and any A&E Al ~(u) we have E(u, O) < E.(u+ A, O). Fix 
x € O, a diffuse hessian D?u € x (0, RY ) such that supp, (D?u(x)) 4 0, a point 
X, € supp, (D?u(x)) and € € R. Recalling (2.3.2), we define the affine map 


A(z) := €@ Fl (x, u(x), Du(x),X,)-(z-2), zeR". 


Fix ¢ > 0 small, x € 2 and choose as O the domain O-.(x) of Lemma 2.2.4.1b). 
Then, tA € An. (wu) for any t € R. Consequently, in view our the above we have 


Biss (u, O.(2)) < Es (u +tA, O.(z)) 


and by applying Lemma 2.2.4.1a), this yields 


0 < max { Hp(z,u(z), Du(2)) : DA(z) + H,(z, u(z), Du(z)) A(z) 


z€O-<(2) 


e—>0 


——> Hp(z, u(x), Du(z)) : (caF (x, u(z), u(x), X.)), 


Hence, 
&- (Hp(e, u(2), Du(x)) Fil (x, u(x), Du(e), X.) J =O, 


for any € € R. By the arbitrariness of € we infer that 


Hp(a, u(x), Du(«)) Fi (x, u(x), Du(x),X,) =0 


for any D?u € Y(Q, Re), a €Q and X, € supp,(D?u(z)), as desired. 


Conversely, let us fix O € 0,2 € O(u), Deu € ¥(Q,RN”), X,, € supp,(D?u(zx)) 
and € € R™ corresponding to a map A € Als (u), Let r be the function of Lemma 
2.2.4.3. By applying Lemma 2.2.4.3 to the above setting, we have 


Dr(0*) > max { Hp(y,u(y),Du(y)) : DA(y) + Hy(y.u(y),Du(y))- AW) 
> Hp(2,u(z),Du(2)) : DA(z) + H,(x,u(a),Du(2)) - A(z) 
= Hp(2,u(2),Du(a)) : (€ @ Fl, (2, u(a),D ks, 
=e (Hp(e,u (x), Du(x)) Fl (x, u(x), Du(x), X.) 


and hence Dr(0*) > 0 because u is a D-solution. Due to the fact that r(0) = 0 
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and r is convex, by inequality (2.2.6) we have r(t) > 0 for all t > 0. Therefore, 


Eoo(t, O) < Ex(u+A,O), VOEQR, V AE Al™(u), 


The case of A € As is completely analogous. Fix D?u € Y (2, RN’), O€EQ, 
a € O(u), X, € supp,(D?u(x)) and an A with A(x) 1 R(Hp(2, u(x), Du(x))) and 
DAE 2a; A(x), X,). By applying Lemma 2.2.4.3 again, we have 


Dr(0*) > max {Hp(y,u(y),Du(y)) : DA(y) + H,(y.u(y),Du(y))- AW) 
> Hp(z,u(z),Du(x)) : DA(z) + H, (x, u(x), Du(x)) - A(z). 


If Hp(x, u(x), Du(x)) 4 0, then by the definition of f(x, A(x), X,) we have 
Dr(0*) > Hp(a, u(x), Du(z)) : DA(x) + H, (x, u(x), Du(zx)) - A(a 
— A(z) - (FS (#,u(x),Du(z),X.) — Hy(2, ae - i) 
= — A(e)"[Hp(x, u(x), Du(e)) + (Fe (#, u(x), Du(w), Xe) 
— H,(x,u(x), Du(2)) ) 


and hence Dr(0*) > 0 because wu is a D-solution on Q. If Hp(x, u(x), Du(x)) = 0, 
then again Dr(0T) > 0 because A(x) = 0. In either cases, by inequality (2.2.6) we 
obtain r(t) > 0 for all t > 0 and hence 


E(u, O) < E,(u+A,O), VOER, VAE AG™(u). 


The theorem has been established. 


2.3.5 Proof of Corollary 2.1.1 


If u € C?(Q, RY), then by Lemma 2.2.2 any diffuse hessian of u satisfies D?u(x) = 
Op2u(z) for ae. © € Q. By Remark 2.3.3, we may assume this happens for all 
« € Q. Therefore, the reduced support of D?u(x) is the singleton set {dp2.(2)}- 
Hence, for Al (u), we have that any possible affine map A satisfies DA = 
D(€H(z, u(x), Du(z))) and A(x) = 0. In the case of As (u), we have that any 
possible affine map A satisfies 


A(x)'Hp(a,u(x),Du(z)) =0, DAE £ (2, A(x), D*u(z)), 
which gives 
DA(z) : Hp(z, u(x), Du(x)) = —A(z) - (Hpp(2, u(x), Du(x)) : D?u(x) + 
) 


+ Hp, (2, u(x), Du(a)) : Du(w) + Hpe(x,u(x), Du(x)) :1) 
= —A(x) - Div(Hp(-, u, Du)) (2). 
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As a consequence, the divergence Div (A'H P(, ah Du)) (x) vanishes because 
DA(x) : Hp(2, u(x), Du(x)) + A(x) - Div(Hp(-, u, Du))(z) = 0. 


The corollary has been established. 
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Chapter 3 


Rigidity and flatness of the image 
of certain classes of mappings 
having tangential Laplacian 


3.1 Introduction 


Suppose that n,N are integers and (2 an open subset of R”. In this paper we 
study geometric aspects of the image u(Q) C R% of certain classes of C? vectorial 
solutions u : R” D Q —+ R% to the following nonlinear degenerate elliptic PDE 
system: 


[DuJtAu=0 in Q. (3.1.1) 
Here, for the map u with components (wy, ..., uy)! the notation Du symbolises the 


gradient matrix 


Du(z) = (Ditte(x)) 7" eR**™ , Dy= 0/er,, 


elm 


Au stands for the Laplacian 


Mie) = S° Dju() € RY 
i=1 


and for any X € R‘*", [X]-+ denotes the orthogonal projection on the orthogonal 
complement of the range of linear map X : R" —> RY: 


[X]}- = Projpyxys- 


Our general notation will be either self-explanatory, or otherwise standard as e.g. 
in [32, 38]. Note that, since the rank is a discontinuous function, the map [-]+ 
is discontinuous on R**”; therefore, the PDE system (3.1.1) has discontinuous 
coefficients. The geometric meaning of (3.1.1) is that the Laplacian vector field 
Au is tangential to the image u(Q) and hence (3.1.1) is equivalent to the next 
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statement: there exists a vector field 


A:R*D0—->R* 


such that 
Au=DuA inQ. 


As we show later, the vector field is generally discontinuous (Lemma 3.2.1). 


Our interest in (3.1.1) stems from the fact that it is a constituent component 
of the p-Laplace PDE system for all p € [2,00]. Further, contrary perhaps to 
appearances, (3.1.1) is in itself a variational PDE system but in a non-obvious 
way. Deferring temporarily the specifics of how exactly (3.1.1) arises and what is 
the variational principle associated with it, let us recall that, for p € [2,00), the 
celebrated p-Laplacian is the divergence system 


A,u := Div([Dul/’-7Du) =0 ind (3.1.2) 


and comprises the Euler-Lagrange equation which describes extrema of the model 
p-Dirichlet integral functional 


E(u) = f (Dur. u € WO, RY), (3.1.3) 


in conventional vectorial Calculus of Variations. Above and subsequently, for any 
X € R**", the notation |X| symbolises its Euclidean (Frobenius) norm: 


N ‘ 1/2 
IX| = pas 2 | 


a=1 i=1 


The pair (3.1.2)-(3.1.3) is of paramount important in applications and has been 
studied exhaustively. The extremal case of p > 00 in (3.1.2)-(3.1.3) is much more 
modern and intriguing, in that totally new phenomena arise which are not present 
in the scalar case. It turns out that one then obtains the following nondivergence 
PDE system 


Agt = (Du @ Du + |Du|?[Du]*@ 1) ‘De=0: im Q, (3.14) 
which is known as the oo-Laplacian. In index from, (3.1.4) reads 
N n 
SS > (Dita Djug + |Du|?[Du]z 5:3 ) Dsus = Dy ty =A yecey NV. 
B=1 i,j=1 


The system (3.1.4) plays the role of the Euler-Lagrange equation and arises in 
connexion with variational problems for the supremal functional 


Ex(u, O) := ||Dullr~@, uu E WY?(O,R), O EQ. (3.1.5) 
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The scalar case of N = 1 in (3.1.4)-(3.1.5) was pioneered by G. Aronsson in 
the 1960s [4-8] who initiated the field of Calculus of Variations in L®, namely 
the study of supremal functionals and of their associated equations describing 
critical points. Since then, the field has developed tremendously and there is an 
extensive relevant literature (see e.g. [16-19, 21, 26, 46, 72, 73] and the lecture 
notes [15, 28, 54]). In particular, although vectorial supremal functionals began 
to be explored early enough, the oo-Laplace system (3.1.4) which describes the 
necessary critical conditions in L°° in the vectorial case N > 2 first arose in the 
early 2010s in [49]. The area is now developing very rapidly due to both the 
mathematical significance as well as the importance for applications in several 
areas (see [2, 14, 31, 36, 63], [50, 52, 53, 56—-58)). 


In this paper we focus on the C? case and establish the geometric rigidity and 
flatness of the images of solutions u : R” D Q —+ R% to the nonlinear system 
(3.1.1), under the assumption that either Du has rank at most 1, or that n = 2 
and wu has an additively separated form, see (3.1.6). As a consequence, we obtain 
corresponding flatness results for the images of solutions to (3.1.2) and (3.1.4). 
Both aforementioned classes of solutions furnish particular examples which provide 
substantial intuition for the behaviour of general extremal maps in Calculus of 
Variations in L®, see e.g. [9, 10, 28, 50, 53, 54, 63] where solutions of this form 
have been studied. Obtaining further information for the still largely mysterious 
behaviour of oo-Harmonic maps is perhaps the greatest driving force to isolate and 
study the particular nonlinear system (3.1.1). For example, it is not yet know to 
what extend the possible discontinuities of the coefficients relates to the failure of 
absolute minimality. 


It is also worth clarifying that, although as it is well-known the Dirichlet prob- 
lem over a bounded domain may not in general be solvable for the oo-Laplacian 
not even in the scalar-valued case, if one does not prescribe boundary values (and 
we do not in this paper) it can be demonstrated that infinitely many non-trivial 
classical solutions do exist, in particular of the form arising in this paper (see for 
instance the explicit constructions of C? solutions in [50]). Therefore, the results 
herein are non-void and numerous solutions as those exhibited herein do exist. 


Let us note that the rank-one case includes the scalar and the one-dimensional 
case (i.e. when min{n, NV} = 1), although in the case of N = 1 (in which the 
single co-Laplacian reduces to Du @ Du : D?u = 0) (3.1.1) has no bearing since it 
vanishes identically at any non-critical point. 


The effect of (3.1.1) to the flatness of the image can be seen through the L® 
variational principle introduced in [52], wherein it was shown that solutions to 
(3.1.1) of constant rank can be characterised as those having minimal area with 
respect to (3.1.3)-(3.1.5). More precisely, therein the following result was proved: 
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3.1.1 Theorem [cf. [52, Theorem 2.7, Lemma 2.2]] 


Given N > n> 1, let wu: R” DQ —> R% be a C? immersion defined on the open 
set 9 (more generally wu can be a map with constant rank of its gradient on 2). 
Then, the following statements are equivalent: 


1. The map uw solves the PDE system (3.1.1) on Q. 


2. For all p € [2,00], for all compactly supported domains O € 2 and all C! 
vector fields vy : O —+ R% which are normal to the image u(O) C RY 
(without requiring to vanish on 0O), namely those for which v = [Du]*v in 
O, we have 


|Dul|z2(0) < ||Du + Dy||z-). 


3. The same statement as in item (2) holds, but only for some p € [2, ov]. 


If in addition p < oo in (2)-(3), then we may further restrict the class of normal 
vector fields to those satisfying v|ao = 0 (see Figure 1). 


In the paper [52], it was also shown that in the conformal class, (3.1.1) expresses 
the vanishing of the mean curvature vector of u(Q). 


The effect of (3.1.1) to the flatness of the image can be easily seen in the case 
of n=1< WN as follows: since 


[w]-u” =0 nQCR 
and in one dimension we have 
u' @u! 
[u’]+ = r= Tue” on {u' # Of, 
i on {uw =O}, 


we therefore infer that wu” = fu’ on the open set {u’ 4 0} C R for some function 
f, readily yielding after an integration that u(Q) is necessarily contained in a 
piecewise polygonal line of RY. As a generalisation of this fact, our first main 
result herein is the following: 


Q 


Figure 1. Illustration of the variational principle characterising (3.1.1). 
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3.1.2 Theorem [Rigidity and flatness of rank-one maps with 
tangential Laplacian] 


Let Q C R” be an open set and n, N > 1. Let u € C?(O,R*%) be a solution to the 
nonlinear system (3.1.1) in Q, satisfying that the rank of its gradient matrix is at 
most one: 


rk(Du) <1 inQ. 


Then, its image u(Q) is contained in a polygonal line in RY, consisting of an at most 
countable union of affine straight line segments (possibly with self-intersections). 


Let us note that the rank-one assumption for Du is equivalent to the existence of 
two vector fields € : R” D Q —> R® and a: R”® D Q —> R" such that Du = €@a 
in Q. 

Example 3.1.3 below shows that Theorem 3.1.2 is optimal and in general rank- 
one solutions to the system (3.1.1) can not have affine image but only piecewise 
affine. 


3.1.3 Example 


Consider the C? rank-one map u : R? —> R? given by 


_ (ee) x < 0, YE R, 
ula, y) = { (t24,c*), c>0, yER 


Then, u=vof with v: R — R? given by v(t) = (t, |t|) and f : R? — R given 
by f(x,y) = sgn(x)x* (see Figure 2). 


Figure 2. The graph of the function f and the image of the curve v comprising u. 


It follows that u solves (3.1.1) on R?: indeed, Au is a non-vanishing vector field on 
{x # 0}, being tangential to the image thereon since it is parallel to the derivative 
v(t) = (1,41) for ¢ 4 0. On the other hand, on {x = 0} we have that Au = 0. 
However, the image u(R?) of u is piecewise affine but not affine and equals v(R). 
Note that (3.1.1) is undetermined, especially without the requirement of boundary 
conditions. Therefore, the point of this example is to show that the solutions in 
general do not have affine image, although some of them may do, for instance the 


30 


trivial affine ones. 


As a consequence of Theorem 3.1.2, we obtain the next result regarding the 
rigidity of p-Harmonic maps for p € [2,00) which complements one of the results 
in the paper [53] wherein the case p = oo was considered. 


3.1.4 Corollary [Rigidity of p-Harmonic maps, cf. [53]] 


Let © C R” be an open set and n,N > 1. Let u € C?(Q, RY) be a p-Harmonic 
map in 2 for some p € [2,00), that is u solves (3.1.2). Suppose that the rank of 
its gradient matrix is at most one: 


rk(Du) <1 inQ. 


Then, the same result as in Theorem 3.1.2 is true. 


In addition, there exists a partition of Q to at most countably many Borel sets, 
where each set of the partition is a non-empty open set with a (perhaps empty) 
boundary portion, such that, on each of these, u can be represented as 


CSI fF, 


Here, f is a scalar C? p-Harmonic function (for the respective p € [2, 00)), defined 
on an open neighbourhood of the Borel set, whilst v : R —> R% is a Lipschitz 
curve which is twice differentiable and with unit speed on the image of f. 


Now we move on to discuss our second main result which concerns the rigidity of 
solutions u : R? D Q —> RN to (3.1.1) for N > 2, having the additively separated 
form 


u(x,y) = f(x) — fly) (3.1.6) 


for some curve f : R —> R%. Solutions of this form are very important in relation 
to the oo-Laplacian. If N = 1, all co-Harmonic functions of this form after a 
normalisation reduce to the so-called Aronsson solution on R? 


_ |x|4/3 ee 


u(x, y) = |y 


which is the standard explicit example of a non-C? oo-Harmonic function with 
conjectured optimal regularity. In the vectorial case, the family of separated so- 
lutions is quite large. For N = 2, a large class of such vectorial solutions was 
constructed in [50] and is given by 


uUsy) = [ (costco). sin( (0) dt 


with K a function in C1(R) satisfying certain general conditions. The simplest 
non-trivial example of an oo-Harmonic map with this form (defined on the strip 
{|x — y| < 7/4} C R?) is given by the choice K(t) = t. Our second main result 
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asserts that solutions of separated form to (3.1.1) have images which are piecewise 
affine, contained in a union of intersecting planes of R“. More precisely, we have: 


3.1.5 Theorem [Rigidity and flatness of maps with tangen- 
tial Laplacian in separated form] 


Let Q C R? be an open set and let also N > 2. Let u : R" D 2 —> R® bea 
classical solution to the nonlinear system (3.1.1) in Q, having the separated form 
u(x, y) = f(x) — f(y), for some curve f € (W?? 7 C?)(R, RY) and some p > 1. 


Then, the image u(Q) of the solution is contained in an at most countable union 
of affine planes in R’. 


In addition, the proof of Theorem 3.1.5 shows that every connected component 
of the set {rk(Du) = 2} is contained entirely in an affine plane and every connected 
component of the set {rk(Du) < 1} is contained entirely in an affine line. 


Note that our result is trivial in the case that N = n = 2 since the codimension 
N —n vanishes. Additionally, due to the regularity of the solutions, if a C? 
mapping has piecewise affine image, then second derivatives must vanish when 
first derivatives vanish at the “breaking points”. Further, one might also restrict 
their attention to domains of rectangular shape, since any map with separated form 
can be automatically extended to the smallest rectangle containing the domain. 


Also, herein we consider only the illustrative case of n = 2 < N and do not 
discuss more general situations, since numerical evidence obtained in [63] suggests 
that Theorem 3.1.5 does not hold in general for solutions in non-separated form. 


In this paper we try to keep the exposition as simple as possible and therefore 
we refrain from discussing generalised solutions to (3.1.1) and (3.1.4) (or (3.1.2)). 
We confine ourselves to merely mentioning that in the scalar case, oo-Harmonic 
functions are understood in the viscosity sense of Crandall-Ishii-Lions (see e.g. 
[28, 54]), whilst in the vectorial case a new candidate theory for systems has been 
proposed in [58] which has already borne significant fruit in [14, 31, 56-58, 63). 


We now expound on how exactly the nonlinear system (3.1.1) arises from (3.1.2) 
and (3.1.4). By expanding the derivatives in (3.1.2) and normalising, we arrive at 


|Dul? 


D 
Du ® Du: D?u + —— Au =i: (3.1.7) 
p— 


For any X € R‘*", let |X]! denote the orthogonal projection on the range of the 
linear map X : R” —> R*: 


[X]! := Projpcx): (3.1.8) 


Since the identity of R% splits as I = [Du]!+[Du]+, by expanding Au with respect 
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to these projections, 


Du Dul? 
Deepa Dag Per plas 
p—2 Pp 


The mutual perpendicularity of the vector fields of the left and right hand side 
leads via a renormalisation argument (see e.g. [49, 52, 53]) to the equivalence of 
the p-Laplacian with the pair of systems 


Du ® Du: D?u + Ba ue [Du]'Au=0, |Du|?[Du]+Au = 0. (3.1.9) 
p— 


The oo-Laplacian corresponds to the limiting case of (3.1.9) as p — oo, which 
takes the form 
Du@®Du:D?u=0, |Du|?[Du]+Au = 0. (3.1.10) 


Hence, the oo-Laplacian (3.1.4) actually consists of the two independent systems 
n (3.1.10) above. The system |Du|?[Du]+Au = 0 is, at least on {Du 4 0}, 
equivalent to (3.1.1). Note that in the scalar case of N = 1 as well as in the case 
of submersion solutions (for N <n), the second system trivialises. 


We conclude the introduction with a geometric interpretation of the nonlinear 
system (3.1.1), which can be expressed in a more geometric language as follows:! 
Suppose that u(Q) is a C? manifold and let A(u) denote its second fundamental 
form. Then 

[Du}* Au = —tr A(u)(Du, Du). 


The tangential part [Du]! Au of the Laplacian is commonly called the tension field 
in the theory of Harmonic maps and is symbolised by 7(w) (see e.g. [69]). Hence, 
we have the orthogonal decomposition 


Au = t(u) — tr A(u)(Du, Du). 


Therefore, in the case of higher regularity of the image of u, we obtain that the 
nonlinear system 


Aa te) We, (3.1.11) 


is a further geometric reformulation of our PDE system (3.1.1). 


3.2 Proofs 


In this section we prove the results of this paper. Before delving into that, we 
present a result of independent interest in which we represent explicitly the vector 
field A arising in the parametric system Au = DuA, in the illustrative case of 
n= 2. 


We will be using the symbolisations “cof” ,“det” and “rk” to denote the cofactor 


'This fact has been brought to our attention by Roger Moser. 
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matrix, the determinant function and the rank of a matrix, respectively. 


3.2.1 Lemma [Representation of A] 


Let u € C?(Q,R*) be given, 9 C R? open, N > 2. The following are equivalent: 
1. The map uw is a solution to the PDE system (3.1.1). 


2. There exists a vector field A : R? D Q —> R® such that 


Au=DuA in. 


In (2), as A one might choose 


cof(Du™Du)" + = 
det(DuTDay D4) "Au on {rk(Du) = 2}, 
Veer ne on {rk(Du) = 1}, 


|Du Du! |? a 
0, on {rk(Du) = 0}. 


A is uniquely determined on {rk(Du) = 2} but not on {rk(Du) < 2} and any 
other A has the form A+ V, where V(z) lies in the nullspace of Du(x), « € Q. 


3.2.2 Proof of Lemma 3.2.1 


The equivalence between (1)-(2) is immediate, therefore it suffices to show that 
A satisfies Au = Du A and is unique on {rk(Du) = 2}. Let A be as in (2). On 
{rk(Du) = 2}, the 2 x 2 matrix-valued map Du' Du is invertible and 


cof (Du™Du)' 


Du! D: ae >, 
( m u) det (Du' Du) 


Since Du'Au = Du' Du A, we obtain that A = A. 


The claim being obvious for {rk(Du) = 0} = {Du = 0}, it suffices to consider 
only the set {rk(Du) = 1} in order to conclude. Thereon we have that Du can be 
written as 


Di=E@a, in{rk(Du)= 1}, 


for some non-vanishing vector fields € and a. By replacing € with €|a| and a with 
a/\a|, we may assume |a| = 1 throughout {rk(Du) = 1}. If Au = DuA, we 
have Au = (€ @ a)A and since any component of A which is orthogonal to a is 
annihilated, we may replace A by Aa for some function A. Therefore, 


Au= (€@a)A = (€@a)(Aa) = EAlal? = E 
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and hence €- Au = A\€|? and also €'Du = al€|?. On the other hand, since 
DuDu' = (€@a)a@gp=E@E, [DuDu"| = EI’ 


we infer that 


Au-€ ar) Au! (€ @ €) Du + DuDu! 
" ( 2 )( | € EP coe 


as claimed. 


We now continue with the proof of the main results. 


The main analytical tool needed in the proof of Theorem 3.1.2 is the next 
rigidity theorem for maps whose gradient has rank at most one. It was established 
in [53] and we recall it below for the convenience of the reader and only in the case 
needed in this paper. 


3.2.3 Theorem [Rigidity of Rank-One maps, cf. [53]]| 


Suppose Q C R” is an open set and u is in C?(Q,R%). Then, the following are 
equivalent: 


(i) The map u satisfies that rk(Du) < 1 on Q. Equivalently, there exist vector fields 
€:0Q —>+R® anda: 2 — R" with a € C'(Q,R”) and € € C1(Q \ {a = 0}, RY) 
such that 


Du=€@a, on?. 


(ii) There exists Borel subset {B;}icn of Q such that 


i=1 


and each B; equals a non-empty connected open set with a (possibly empty) bound- 
ary portion, functions {f;}ien € C?(Q) and curves {v;}ien C Wig (R, RY) such 
that, on each B; the map u has the form 


U=, 6 Ff; on Be (3.2.1) 


Moreover, |v;| = 1 on the interval f;(B;), v; = 0 on R \ f,(B;) and vy’ exists 
everywhere on f;(B;), interpreted as 1-sided derivative on Of;(B;) (if f;(B;) is not 
open). Also, 


Du = (Yo fi) @Dfi , on B;, (3.2.2) 


Du = (Yio f) @DA@Df+ of) @D' fi, on Bi. 
In addition, the local functions (f;)% extend to a global function f € C?(Q) with 
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the same properties as above if is contractible (namely, homotopically equivalent 
to a point). 


We may now prove our first main result. 


3.2.4 Proof of Theorem 3.1.2 


Suppose that wu: R” D Q —> R% is a solution to the nonlinear system (3.1.1) in 
C?(Q,R%) which in addition satisfies that rk(Du) < 1 in Q. Since {Du = 0} is 
closed, necessarily its complement in Q which is {rk(Du) = 1} is open. 


By invoking Theorem 3.2.3, we have that there exists a partition of the open 
subset {rk(Du) = 1} to countably many Borel sets (B;)$° with respective functions 
(fi)f° and curves (v;)?° as in the statement such that (3.2.1)-(3.2.2) hold true and 
in addition 

Df; 4 0 on B;, 1 EN. 


Consequently, on each B; we have 


(Yi0 fi) @ (42 fi) 


\Y%° Fil? 


[Du]+ = [Mo f,) @DfJ+ =1- 


Au = (fo fi)|Dfi? + Ho fo Afi. 


Hence, (3.1.1) becomes 


(vio fi) ®@(Y%0 fi) 


_ 
Yio fil? 


(ere mIDAP + Wo MAK) =o, 


on B;. Since |v;|? = 1 on f;(B;), we have that 1’ is orthogonal to v/ thereon and 
therefore the above equation reduces to 


(v;' 0 f)|DFiP = 0 onB, ieN. 


Therefore, v; is affine on the interval f;(B;) C R and as a result u(B;) = v;(f;(B;)) 
is contained in an affine line of R’, for each i € N. On the other hand, since 


u(Q) = u({Du = OF) LJ u(Bi) 


iEN 


and u is constant on each connected component of the interior of {Du = 0}, the 
conclusion ensues by the regularity of u because u({Du = O}) is also contained in 
the previous union of affine lines. The result ensues. 


Now we establish Corollary 3.1.4 by following similar lines to those of the re- 
spective result in [53]. 


3.2.5 Proof of Corollary 3.1.4 
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Suppose uw is as in the statement of the corollary. By Theorem 3.2.3, there exists, 
a partition of 2 to Borel sets {B;}ien, functions f; € C?(Q) and Lipschitz curves 
{vi}icn : R — R® with |v{] = 1 on f,(B;), || = 0 on R \ f;(B;) and twice 
differentiable on f;(B;), such that ulp, = 1% 0 f; and (3.2.2) holds as well. Since on 
each B; we have 


[Du| = |(%° fi) ® DF] = [Dfil, 
by (3.1.7) and the above, we obtain 


(vic fi) @Dfi) ® (Ho fi) @Df) : tu o f;) @Df, @Df; + (Vio fi) ® D’ f, 


ID fil? 


i =) 


c of) Aji a (© pips} =. 0; 
on B;. Since vj’ is orthogonal to 1; and also v/ has unit length, the above reduces 
to 
! ~T2 ID fil? i " 4 _ 
(¥,0 fi)|Df,® Df, : D°fi + —S AR] + —5z% eo f)IDAI" = 9, 
p—2 p-—2 
on B;. Again by orthogonality, the above is equivalent to the pair of independent 
systems 


Df;|? , 
10 f) [Phe Df Diet oon =0, (Wof)|DAily = 0, 
on B;. Since |v{| = 1 of f;(B;), it follows that A,f; = 0 on B; and since (B;)?° 
is a partition of 2 of the form described in the statement, the result ensues by 


invoking Theorem 3.1.2. 


We may now prove our second main result. 


3.2.6 Proof of Theorem 3.1.5 


The system [Du]+Au = 0 is equivalent to 
Au =Du-A (3.2.3) 
for a vector field A with components a,b. Then (3.2.3) can be rewritten as 
f'(e) — Fy) = ala, y) f(x) — o(@, y) Fy). (3.2.4) 
The choices (x, y) = (z,z +t) and (x,y) = (z +t, z) in (3.2.4) yield the equations 
f"(2) —f"@t+t) = a(z,z+bf'(z) — b(z,24+t)f'(z +t) (3.2.5) 


and 
f"(<4+t) -—f"(2%) = az +t, 2) f'(e +t) — b(z +t, z) f(z) (3.2.6) 


37 


respectively. Let f, denote the a-component of f, a = 1,...,N.. By subtracting 
(3.2.5) from (3.2.6) we get for t 4 0 that 


9 Falz +t) ~ fal) falz +t) — fal) 
t 


= (a(z +t,z) + b(z,2+ i)) ; 


F a(z+t,z)—a(z,z+t)  bd(z,2+t)-—d(z+t,z) 
+ fl2)( ; + ; ) 


(3.2.7) 


for a =1,..., N. On the set {f% = 0}, equation (3.2.7) becomes 


2fn (2) = (@(z,2) + 6(2,2)) fa (2) (3.2.8) 


as t + 0. Note also that {f' = 0} is closed and its complement { f’, 4 0} is open. 
Now let us set 


Cy(z,t) = — e esta ls) 


On {f, 4 0}, (3.2.7) yields that 
falz +t) — fa(z) 
t 


C.40 2 ae E 


Fix an index a € {1,..., N}, 6 > 0, an infinitesimal sequence (t,,)?° and consider 
the inner 6-neighbourhood O; of the set {f/, # 0}, namely 


_ (a(ett, 2) +d(z, z+) fale) f) 


t 


Oy = {« ER”: fi(x) #0 and dist(x, Off, 4 0}) > 5}. 


Then for any fixed 6 > 0 small, there exists a constant cs > 0 such that along the 
sequence t,, — 0 we have 


lL falta) = Je0) 
ta) Co 


[Cab ta)lvon £2 


LP(O5) 
Ll dain) = ie) 
+ lla + b| L&°(Q) fi o : 
a m LP(O5) 


(3.2.9) 


IA 


1 
a (a1 |z7(05) + |la+ bul lie 


IA 


1 
= (217 |z>(py + |la + bulla 


Note that the right hand side of the above estimate is bounded uniformly in m € N 
as f’” € L?(R,R*) and f’ € C'(R,RY). By letting 6 > 0 and using a standard 
diagonal argument, (3.2.9) implies that there exists a function C, such that 


Cal=sten) — Co in de # O}), 
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as m —> oo along a subsequence of indices (m;)¢°. As a result, (3.2.7) becomes 


2f0 (2) = (Gz, 2) + W(z,2)) falz) + f(2)Ca(z) on {f, #0}, (3.2.10) 


for any a = 1,..., N. Combining equations (3.2.8) and (3.2.10), we infer that there 
exist measurable functions A, B : R —> R such that 


f" =Af'+ Bf" ae. onR. (3.2.11) 


The goal in now to show that (3.2.11) implies that the torsion of the curve f 
vanishes, at least on a union of subintervals of IR. The idea to project on three- 
dimensional subspaces of R% in order to utilise standard ideas of elementary dif- 
ferential geometry of curves. 


To this end, let P; : RY —+ R% be the orthogonal projection on a 3D subspace 
V3 = P3(R%) of R%. The choice of 3-dimensional subspaces owes to the fact that 
we would like to use the classical formulas of differential geometry of curves in the 
Euclidean space. Then, P3f : R —> V3 & R? is a curve in R°, which is C?. By 
(3.2.11) we have, 


(P3f)" = A(P3f)’ + B(P3f)” a.e. on R. 


Let "x" denote the cross (exterior) product in R°. Then, the curvature of P3f is 
given by 
K = |(P3f)' x (P3f)"| 
and, on {k 4 0}, the torsion is given by 
[(P3f)! x (Ps if)" (Ps fy" 
\(Psf)! x (Pf? 


Note that {« 4 0} is open, as P3f is C?. Then, we have: 


e On the topological interior int({« = 0}), P3f is contained in an affine line of V3. 


e On the topological interior int({« 4 0}), P3f is planar and hence contained in 
affine plane of V3. 


Since O({« = 0}) is nowhere dense, it follows that f(O({« = 0})) is contained 
in the boundary of an affine plane or an affine line. Hence, we have that, for any 
projection P3f on a 3D subspace of R%, the projected curve is contained in an at 
most countable union of affine planes and lines. Therefore, the same is true for 
f itself by elementary analytic geometry: if all 3-dimensional projections of the 
image set in the space R% for n > 3 are planes or lines, the same is true for the 
image itself. The conclusion follows. 
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Chapter 4 


Explicit oo-harmonic functions in 
high dimensions 


4.1 Introduction 


Let Q C R” be an open set and u € O7(Q) a continuous twice differentiable 
function. In this paper we study the existence of solutions to the PDE 


Agu = > Du Dju D5 u = 0 (4.1.1) 


ij=l 
of the form 


u(x) = [] Ao, 


where f; are possibly non-linear for 1 < i < n, and x = (q},...,2n)', c € 2. 
Solutions of this form are called separated coo-harmonicfunctions. In the above 
De and D?, = se. The equation (4.1.1) is called oo-Laplacian (being a 
special case of the so-called more general the Aronsson equation) and it arises in 
Calculus of Variations in L® as the analogue of the Euler-Lagrange equation of 
the functional 


E(u, O) := ||Dul|zo), OEQ, we We°(Q,R). 
These objects first arose in the work of G. Aronsson in the 1960s (see [6],[7]) 
and nowadays this is an active field of research for vectorial case N > 2 for u € 
W.2%(Q,.R) which has begun much more recently in 2010s (see e.g. [49]). Since 
then, the field has been developed enormously by N. Katzourakis in the series of 
papers ([50-53, 55-59]) and also in collaboration with the author, Abugirda, Croce, 
Manfredi, Moser, Parini, Pisante and Pryer ((14], [2], [31], [60], [61], [62], [|63-65]). 
A standard difficulty of (4.1.1) is that it is nondivergence form equation and since 
in general smooth solutions do not exist, the definition of generalised solutions is 
an issue. To this end, the theory of viscosity solutions of Crandall-Ishii-Lions is 


AO 


utilised (see e.g. [54]). 


In this paper all the separated coo-harmonicfunctions are found for n = 2 in 
polar coordinates and for all n > 2 in cartesian coordinates. Some of these new 
solutions derived herein coincide with previously known classes of solutions. For 
instance, the well-known G. Aronsson’s solution u(x, y) = |x|3 — |y|? which has a 
C1/3 regularity, described in Remark 4.2.2. Also M.-F. Bidaut-Veron, M. Garcia 
Huidobro and L. Veron have found solutions in [20] which coincide with first two 
solutions of the theorem 4.1.1. In addition I.L. Freire, A. C. Faleiros have found 
solutions of (4.1.1) in [44], but only one of their non-trivial solutions coincides with 
a particular case of Theorem 4.1.2 when A = 1. There may exist other additional 
solutions but this topic is not discussed herein. 


The main results of this paper are contained in the following theorems. 


4.1.1 Theorem [Separated two - dimensional oo - harmonic 
functions in polar coordinates] 


Let u : Q C R? —> R be a C?(Q) separated co-harmonicfunction of the oo- 
Laplace equation in polar coordinates 


2 1 1 
Un Ure + 2 Urtetirg + —j ua too — “gure =) (4.1.2) 


of the form u(r,#) = f(r)g(@). 


(i) Assume |f(r)| = r4 and |g(0)| = e?°, where A and B are any constants, 
then 
A?— A+ B?=0 


or 


(ii) Assume |f(r)| = r4 and |g(0)| = 1g(80) jel*o C4" then G satisfies the 
following 


— arctan oe) + Ara arctan woes if A2-A>O 
1 ‘ 
ae i A= 
t+c= aa : 
— arctan G(t), A= 


GG) A-1 G(t)—VA= A? 
~e aee in| Oe 


, if A?7-A<0O, 


where c is any constant, provided RHS is well defined. 
or 


r &(t 
(iii) Assume |g(@)| = e?? and |f(r)| = |f (ro) let ea where ©® satisfies the 


Al 


following 


1 $?(t)+B? i ot B(t)—5 2 1 
2 In PHS H1B2|— 2 Bt arctan ae if B*—=->0 
1 ©?(t)+B? 1 1 5 2 1 

In |¢| + C= 9 In ©2(t)—® t)+B? 2 B(t)—47 if B AS oe 0 
1 ©? (t)+B? 1 (t)—4—,/4—B? gion 
3 In &2()-B)+B2| 1 1p In B14 IR F if B+ — 4 < 0, 


where c is any constant, provided RHS is well defined. 


4.1.2 Theorem [Separated two-dimensional oo - harmonic 
functions] 


Let u: Q C R? —>+ R be a C?(Q) separated oo-harmonicfunction of the oo- 
Laplace equation 
UsUee + 2UgtyUry + UsUyy = 0 (4.1.3) 


of the form u(x, y) = f(x)g(y). Then, one of the following holds: either 


(i) [f(x)| = |F(ao)|e4@- and |g(y)| = |a(yo) lel CO, where G satisfies 


1 : ~ 
i. OE P . tA= 
ES 1 Git Git : 
—5, arctan =" 4 : (424.620) , otherwise 


or 


(ii) |f(x)| = | f (xo) elo" and |g(y)| = |g(yo)|e24-™ , where F satisfies 


1 : 
FO? if B=0 
Ft) | Ft) 


ti+c= 1 
S7p a 2(B2+F2(t)) 


, otherwise. 


4.1.3 Theorem [Separated n-dimensional oo - harmonic func- 
tions] 


Let n > 2 andu:Q CR" —>R bea C?(Q) separated oo-harmonicfunction of 
the oo-Laplace equation 


"Diu DjuDjju = 0. (4.1.4) 
i,j=l 
Then 
|fi(as)| = |fela?)| AO m*) for 1 <i xj <n 
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and 


” R(t) de 
|fi(as)| = |Fj(29) f ; 


where F; satisfies 


t+c= ee ee arctan Ee) AY : 
(Day (Lay Lae Fe) 
iAj iAj iAj 


4.2 Proofs of main results 


In this section we prove our main results. The general idea of our method, which is 
essentially the same for all our proofs, is to use a substitution to derive a “better” 
PDE. Then, we take any points from the domain which are different only in one 
component put them to the “better” PDE and subtract these two equations from 
each other to get a “new” PDE. 


4.2.1 Proof of Theorem 4.1.1 


We can assume that u 4 0 since if wu is a solution then u+c is also a solution then 
the equation (4.1.2) can be written as 


U2 Urr 2 Up Ug Ure  luzug 1 Up Ue 


_ = 0. 4.2.1 
u2 wu r2uu U r4u2 u r? wu ur ( ) 


Let F = “ and G= “, then F, + F? = “*, Gg + G? = “ and $Fy + $G,4 
FG = “*. Note that u(r,@) = f(r)g(@), hence F does not depend on @, since 


Fir.o) = a Analogously G does not depend on r, since G(r,@) = oe Thus 
(4.2.1) becomes 


2 1 1 1 
F°R, + F°+ SFG? + 3@’Got+ GG* — SFG’ =0. (4.2.2) 
‘a r r ‘i 
Set @ = Fr, then r®, — ® = F,r?. Multiplying (4.2.2) by r*, we have 
(6? + G’)(G? + G? — 6) + rb, + G’G, = 0. (4.2.3) 
We have the following 4 cases for the functions ® and G: 
Case (A) ® and G are constant functions. 
Case ® is constant and G is non-constant functions. 


(B) 
Case (C) ® is non-constant and G is constant functions. 
Case (D) ® and G are non-constant functions. 
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Case (A) Let ® = AandG = B, then (4.2.3) gives A= B = 0 or A?—A+B? = 0 
which can be rewritten as (A—})?+B? = } and as the consequent of substitutions 
f(r) =r4 and g(0) = e® up to a constants. 


Case (B) Let ® = A, then G is a non-constant function satisfying (4.2.3) 


(A? + G?)(A? + G? — A) + G’Gs = 0. (4.2.4) 
Therefore iG 
(A? + G?)(A? + G? — A) = hr 
Consequently 
—G? 
do = d 
/ e | a@s@arec - (4.2.5) 
=f = ac— [ i 7 
— f A+ A2—-A+G?2 
Hence 
— arctan 1) + Ars arctan wrt if A2Z-A>O 
Can if A=0 
t+c= a : 
— arctan G(t), if A=1 
G(t) A-1 G(t)—VA—A? | se 2 
— arctan =* + sfzte in| OEE , if A*-A<0O. 


Case (C) Let G = B, then © is a non-constant function satisfying (4.2.3) 
(©? + B?)(©? + B? — 6) + 86, =0. (4.2.6) 


Therefore re 
(G* + B?)(? + B? — 6) = S70": 
Tr 


Consequently 


fea] ar d® 
r J (62+ B2)(62 — @ + B?) 


= av | aie ao | z dd 
a @2 4 B2 (6—4)2+ B?-+ (®— 4)?+ Bet . 


(4.2.7) 
Hence 

1 ®?(t)+B? 11 o(t)—4 ; 

7 In @2(t)—O(t)+B2| 2 Goa Bra ib? = ag 

il ?(t)+B? 1 1 ' 7 
In|{| +e= 9 In PO) +B? | | 2 B(t)—5? if BY — ga 0 

1 ©? (1) +B? 1 ®(t) —B?| , 

zn lare—sre| — yim |sqcin fim) FB — a <0 


Case (D) Let ® and G be non-constant functions, then there exist r; 4 rp and 
6, A @ such that ®(r,) 4 ®(r2) and G(0,) 4 G(62) satisfying (4.2.3). Thus 


ry ®(171)"®, (71) —®* (11) +B (r1)+28(171)’G(8)"+G(8)"Go(8}+G(0)"— “it 
r2®(r2)"®,(r2)—* (rz) + B(r2)+28(r2)"G(8)"+G(8)"Go(0)+G(8)— ae 
Subtracting (4.2.8) and (4.2.9) we get for any 0 


G?(8)(®(r1) — B(r2))(2(B(r1) + G(r2)) — 1) = re&*(r2),(r2) — B(r2)? + O(r2)* 
— 7,0? (r1)®,(r,) + O(r,)? — B(r,)*. 
(4.2.10) 


Let’s consider two cases. 


Case (I) If there exists r; A rg such that 2(®(r1) + ®(r2)) —1 4 0, then (4.2.10) 
gives that G?(0) is a constant, hence G(@) is a step function. 

Case (II) For any r; # ro we have 2(®(r,) + ®(r2)) — 1 = 0, hence ®(r) is a 
step function. 

For both cases we have a contradiction to C+ regularity for oo-Harmonic func- 
tions in two dimensions (see [39], [75]), since ®(r) = +“ and G(@) = “2 have to 
have at least C°* regularity. 


fi fr 2@) 
Finally integrating £ ' =i ¢ 7 = Gand substituting we get | f(r)| = |f (ro) erro “ 
and |g(0)| = Ig(O)|e!% & (Oat which completes the proof. 0 


4.2.2 Remark [The Arronson solution] 


Let A = = in the Theorem 4.1.1ii, then A? — A > 0 and function G satisfies 


1 3 
— arctan ~G(t), 


3 
t = — arctan —G(t 
+c arctan 7 (i) +5 5 


which can be rewritten as 
27G°(t) + 54G?(t) tan 2(t + c) + 32 tan 2(t + c) =0. 
Solving a third degree equation with respect to G(t), we get 


4tan3(t +c) + tan3(t +c) + tan(t +c) 
3 1 — tan?(¢ + c) 


Ga == 


Therefore 


(1 — tans (¢ + c)) (1+ tans ( cape 3]. 
(tan3(t +c) — tan3 (t+c) +1)? 


[ewa=in 
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Hence 


off, ayar_ [1 — tan3 (0 + o)|]1 + tan3(0 + 0)/3_ 
| tan3(@ +c) — tan3(@ +c) +1)3 
1—tan3(0+¢ 

_ HatatO+9l i, 
|1 + tan*(@ + c)|3 
= | cos3 (0 +ce)- sins (0 + c)| - (8). 


We can ignore c(09) since if c,u + co is a solution then uw is also a solution. 


Finally 


19(9)| = |9(40) |e! Sat 
= |9(9)|(| cos (6 + c) — sin3 (0 +c) 


), 


4.2.3 Remark [The Aronsson solution] 


Let A = —} in the Theorem 4.1.1ii, then A? — A > 0 and function G satisfies 
3 
t+ c = arctan 3G(t) — 2 arctan 5Glt). 


Carrying out a similar series of calculations as in Remark 4.2.2 we can find that 
1 
possible solutions are f(r) = r~3 and g(0) = cos3 (#4) —sin3 (4°), hence u(r, 0) = 


2 2 
O+c O+c 


773 (cos3 (242) - sins (24¢)) is the solution of the oo-Laplace equation which was 


described in [10] as u(r,@) = r~3g(0), where 


t tant 
Hs @=t—2aretan (S2*), ae" 
(1 + 3cos?t)3 2 2 2 


The key fact these two solutions are identically equal is tan § = — tan? 7 


4.2.4 Proof of Theorem 4.1.2 


It is a particular case of the Theorem 4.1.3. 
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4.2.5 Proof of Theorem 4.1.3 


We can assume that u 4 0 since if wu is a solution then u+c is also a solution then 
equation (4.1.4) can be written as 


“. Dyu Dju D?,u 

y> eae = 0. (4.2.11) 
U U U 

ij=l 


Pui Thus (4.2.11) becomes 


U 


Let F; = 2 then D;F; + F? = 2#" and FF; = 


U U 


(> rte) a 3 F}(x;)DiFi(a,) = 0. (4.2.12) 


i=l 
Since u(x) = I] fi(a;), then F; depends only on 2;, consequently 
i=1 
D,F;(2;) = Fj (zi). 


1 2 if. 1 2 2 
Set 2°, 2° € O such that 2° = (41, £2,..., 27, +) En) and 2° = (21, 2a, -.-, £5, ---; En), 


where aj # x* in (4.2.12) and subtract these two equations. We find 
(F#(a}) — FF(@})) (2) FR (es) + 2F P(e} + 2FF@})) + 
tA 
2/1 1 2/2) TM (op 2\ 
Fs (a5) Fj (2;) _ Fy (x5) F (x5) = 0, 


assuming F?(aj) A F?(x7), we have 


297 F(a) = — (x5) ie * (25) (a) — 2F? (a!) — 2F?(x}). (4.2.13) 


iAj j 
LHS of (4.2.13) does not depend on «7; and 2; so 
>) F(a) =e 
1Aj 


for all x;. Then F;(x;) = A;, where A; is a constant for all 1 <i #j <n and 
0 


hence | f;(a;)| = | f;(x°)| e4**-*).. Thus (4.2.12) gives 


(m4 a ee) + Fj (x3) Fj (23) =0, 


iAj 


AT 


consequently 


F? 
(a== Fak, 
(x A? + #3) 
ij 
[, Ba 
hence | f;(x;)| = |f;(#})|e * , where F(t) satisfies 


1 F;(t) 
t+c= ————— arctan il if SA} 40. 


2 [S7 a os (A+R) 
iAj ixj a 


Otherwise (i.e. if SA = 0) 


tAl 
F (a5) + F} (25) Fj(a,) = 0, 
so 
F?(a;) + Fi(aj) = 0, since we assume F?(x;) # F?(x%). (4.2.14) 
Solving (4.2.14) we get Fj(x;) = —.. Hence |fi(xi)| = c; for all i # j and 


wj+c 
|f;(25)| = ¢;(|zj; + c]), where c and c; are constants for alll <i <n. 
If there is no j such that F?(a;) A F?(x7) then F?(a;) = c; for alll <j<n 
and (4.2.12) gives that c; = 0 for all 1 < 7 < n. So |f,(a:)| = Ci, where C; is 
constant for all 7. 


4.3 Numerical approximations of oo - harmonic 
functions 


In this section we illustrate the oo-Harmonic functions derived earlier, depending 
on the parameter(s)(see Figure 4.1 - 4.4). The results illustrate that we may 
have a family of solutions depending on the 27-interval even if the parameter(s) 
is/are fixed. For example: the solution on Figure 4.2h is a combination of those 
in Figure 4.21 and Figure 4.2] when @ belongs to Ist and 2nd 27- interval of the 
domain respectively. Colours are linear colour scaled from minimum to maximum. 
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Figure 4.1: The approximation to u of the Theorem 4.1.1 i, depending 
on the parameters A and B. 


(a) A= 0.25,B = 0.433, (b) A=0.5,B =0.5, A= 0.75, B = 0.433, 
min = 79, max = 1209 min = 66, max = in min = 57, max = 866 
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(d) A = 0.25, B = —0.433, (ec) A=0.5,B =-0.5, — (f) A = 0.75, B = —0.433, 


min = 5, max = 79 min = 2, max = 66 min = 3, max = 57 
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Figure 4.2: The approximation to u of the Theorem 4.1.1 ii, depending 
on the parameter A. 


(a) A= 4/3, (b) A=1.15, (c) A=1, 


min = 0, max = 6 min = 0, max = 7 min = 0, max = 10 


Ga=1s, (ec) A=0.15, (f) A=0, 


min = 0, max = 3 min = 0, max = 1.8 min = 1, max = 39 


oe = 2 2 + 8 


05 


(g) A= —0.15, (h) A = —0.05, 


min = 1, max = 41 min = 1, max = 50 


" | 


as 3S 8 8 & B82 8B 


(i) A = —0.05, (j) A = —0.05, 


min = 1, max = 50 min = 1, max = 48 


05 


50 


Figure 4.3: The approximation to u of the Theorem 4.1.1 iii, depending 
on the parameter B. 
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(a) B= -1/3, (b) B= -1/2, (c) B=0, 
min = 0.12, max = 20 min = 0.04, max = 12 min = 1, max = 40 
(a) B=1/3, (°) B=1/2, - 


min = 1, max = 163.58 min = 1, max = aan min = oa max = 2304 


Figure 4.4: The approximation to u of the Theorem 4.1.2 i, depending 
on the parameter A. 
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(a) A= —0.5, (b) A = —0.25, (c) A = —0.05, 
min = 0.0067, max = 663 min = 0.0821, max = 76 min = 0.0665, max = 21 


(d) A=0, (e) A= 0.05, (f) A =0.25, 
min = 1, max = 100 min = 0.6065, max = 21 min = 0.0821, max = 76 
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Chapter 5 


Vectorial variational principles in 
L® and their characterisation 
through PDE systems 


5.1 Introduction 


Let n,N € N and H € C?2(Q x RY x pe) with Q C R” an open set. In this 
paper we consider the supremal functional 


E.o(u,O) := ess supH(-,u,Du), we Wer(Q;R%), OE Q, (5.1.1) 
O 


defined on maps u : R" D Q —+ R%. In (5.1.1) and subsequently, we see the 
gradient as a matrix map Du = (Dju,)¢5);% : R” DO — R**". Variational 
problems for (5.1.1) have been pioneered by Aronsson in the 1960s in the scalar case 
N = 1 ({4]-[8]). Nowadays the study of such functionals (and of their associated 
PDEs describing critical points) form a fairly well-developed area of vivid interest, 
called Calculus of Variations in L°°. For pedagogical general introductions to the 


theme we refer to [13, 28, 54]. 


One of the main difficulties in the study of (5.1.1) which prevents us from util- 
ising the standard machinery of Calculus of Variations for conventional (integral) 
functionals as e.g. in [34] is that it is non-local, in the sense that a global min- 
imisers u of E(-, 2) in Wea, R”) for some fixed boundary data g may not 
minimise E,,(-,O) in W}°(O;R%). Namely, global minimisers are not generally 
local minimisers, a property which is automatic for integral functionals. The rem- 
edy proposed by Aronsson (adapted) to the vector case is to build locality into the 
minimality notion: 
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5.1.1 Definition [Absolute Minimiser] 


Let u € W°°(0; RY). We say that u is an absolute minimiser of (5.1.1) on Q if 


loc 


¥YOeR, 
VoEW™(O; 


RY) \ = Ex~(u,O) < Exn(ut ¢, 0). (5.1.2) 


In the scalar case of N = 1, Aronsson’s concept of absolute minimisers turns out 
to be the appropriate substitute of mere minimisers. Indeed, absolute minimisers 
possess the desired uniqueness properties subject to boundary conditions and, 
most importantly, the possibility to characterise them through a necessary (and 
sufficient) condition of satisfaction of a certain nonlinear nondivergence second 
order PDE, known as the Aronsson equation ([3, 13, 15-18, 25-27, 29, 48, 68, 78]). 
The latter can be written for functions u € C?(Q) as 


Hp(-,u, Du) - D(H(-, u, Du)) = 0. (5.13) 


The Aronsson equation, being degenerate elliptic and non-divergence when for- 
mally expanded, is typically studied in the framework of viscosity solutions. In 
the above, Hp,H,,H, denotes the derivatives of H(x,n,P) with respect to the 
respective arguments and “.” is the Euclidean inner product. 


In this paper we are interested in characterising appropriately defined minimis- 
ers of (5.1.1) in the general vectorial case of N > 2 through solvability of associated 
PDE systems which generalise the Aronsson equation (5.1.3). As the wording sug- 
gests and we explain below, when N > 2 Aronsson’s notion of Definition 5.1.1 is 
no longer the unique possible L° variational concept. In any case, the extension 
of Aronsson’s equation to the vectorial case reads 


Hp(-, u, Du) D(H(-, u, Du)) 


+ H(-,u, Du) [Hp(-, u, Du)}* (Div (He(-, Du)) — H,(-,1,Du)) ap. (5.1.4) 


In the above, for any linear map A : R" —> R%, [A]+ symbolises the orthogonal 
projection Projp(4): on the orthogonal complement of its range R(A) C RY. We 
will refer to the PDE system (5.1.4) as the “Aronsson system”, in spite of the fact 
it was actually derived by N.Katzourakis in [49], wherein the connections between 
general vectorial variational problems and their associated PDEs were first studied, 
namely those playing the role of Euler-Lagrange equations in L°°. The Aronsson 
system was derived through the well-known method of L?-approximations and is 
being studied quite systematically since its discovery, see e.g. [49]-[50], [57, 66]. 
The additional normal term which is not present in the scalar case imposes an 
extra layer of complexity, as it might be discontinuous even for smooth solutions 
(see [50, 53)). 


For simplicity and in order to illustrate the main ideas in a manner which min- 
imises technical complications, in this paper we restrict our attention exclusively to 
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regular minimisers and solutions. In general, solutions to (5.1.4) are nonsmooth 
and the lack of divergence structure combined with its vectorial nature renders 
its study beyond the reach of viscosity solutions. To this end, the theory of D- 
solutions introduced in [57] and subsequently utilised in several works (see e.g. 
[14, 31, 56, 57]) offers a viable alternative for the study of general locally Lipschitz 
solutions to (5.1.4), and in fact it works far beyond the realm of Calculus of Vari- 
ations in L°°. We therefore leave the generalisation of the results herein to a lower 
regularity setting for future work. 


Additionally to absolute minimisers, for reasons to be explained later, in the 
paper [52] a special case of the next L® variational concept was introduced (therein 
for H(x,9, P) = |PI?): 


5.1.2 Definition [oo-Minimal Map] 


Let u € C1(Q;R%). We say that u is an oo-minimal map for (5.1.1) on Q if (i) 
and (ii) below hold true: 


(i) wis a rank-one absolute minimiser, namely it minimises with respect to essen- 
tially scalar variations vanishing on the boundary along fixed unit directions: 


YVOE, VEER 
V 6 € Ch(O; span{é]) 


\ = E,(u,O) < E.(ut+¢,0). (5.1.5) 


(ii) wu has oo-minimal area, namely it minimises with respect to variations which 
are normal to the range of the matrix field Hp(-,u, Du) and free on the boundary: 


YOEN, VdEC (RR) 


with ¢"Hp(-,u, Du) =0 on O \ = E,(u,O) <E.(ut+¢,O). (5.1.6) 


In the above, 


Co(O;R®) := {ye C'(R"; RY) : P=0 on OO}. 


Note also that when N = 1 absolute minimisers and co-minimal maps coincide, 
at least when {Hp = 0} C {H = 0}. Further, in the event that Hp(-,u, Du) has 
discontinuous rank on O, the only continuous normal vector fields ¢ may be only 
those vanishing on the set of discontinuities. 


In [52] it was proved that C? oo-minimal maps of full rank (namely immer- 
sions or submersions) are oo-Harmonic, that is solutions to the so-called oo- 
Laplace system. The latter is a special case of (5.1.4), corresponding to the choice 
H(z,7, P) = ake 


Du D(|Du|?) + |Du|? [Du]*Au = 0. (5.1.7) 
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The fullness of rank was assumed because of the possible discontinuity of the 
coefficient [Du]+, which may well happen even for smooth solutions (for explicit 
examples see [50]). In this paper we bypass this difficulty by replacing the orthog- 
onal projection [-]+ by the projection on the subspace of those normal vectors 
which have local normal C extensions in a open neighbourhood: 


5.1.3 Definition [Orthogonal Projection] 


Let V: R° DQ — R**” be a matrix field and note that 


where for any x € 9, N(V(x)') is the nullspace of the transpose V(r)! € R"*%, 
We define the orthogonal projection 


[V(z)}° = Projgy(eyt)s [V)}-: R7DQ—>R™*, 


where N(V(z)') is the reduced nullspace, given by 


N(V(2)") := {é EN(V(z)") |Je>0 & JECC'UR™ RY): 


E(x) = € & Ely) EN(V(y)"), Vy € B.(a)}. 


It is a triviality to check that N(V(x)") is indeed a vector space and that 


where [V (x)|* = Projyyy(q)t). Note that the definition could be written in a more 
concise manner by using the algebraic language of sheaves and germs, but we 
refrained from doing so as there is no real benefit in this simple case. 


The first main result in this paper is the next variational characterisation of 
the Aronsson system (5.1.4). 


5.1.4 Theorem [Variational Structure of Aronsson’s sys- 
tem] 


Let u: R" D OQ —> R® be a map in C?(Q; RY). Then: 


(1) If uw is a rank-one absolute minimiser for (5.1.1) on Q (Definition 5.1.2(i)), then 
it solves 
Hp(-,u, Du) D(H(-,u, Du)) = 0 on Q. (5.1.8) 


The converse statement is true if in addition H does not depend on 7 € R™ and 
Hp(-, Du) has full rank on 2. 
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(II) If w has co-minimal area for (5.1.1) on Q (Definition 5.1.2(ii)), then it solves 
H(-, u, Du) [Hp(-,u, Du)}+ (Div (Hp(-,, Du)) — Hy(-,u, Du) ) = 0 on. (6.1.9) 


The converse statement is true if in addition for any x € 2, H(z,-,-) is convex on 
R” x RN xn, 


(III) If uw is oo-minimal map for (5.1.1) on Q, then it solves the (reduced) Aronsson 
system 


Axou := Hp(-,u, Du) D(H(., u, Du)) 
+ H(.,u, Du) [Hp(-, u, Du)}+ (Div (Hp(-,u, Du)) = H,,(-,u,Du)) = 


The converse statement is true if in addition H does not depend on 7 € R¥, 
Hp(-, Du) has full rank on Q and for any x € 2 H(z,-) is convex in RN*". 


The emergence of two distinct sets of variations and a pair of separate PDE 
systems comprising (5.1.4) might seem at first glance mysterious. However, it is 
a manifestation of the fact that the (reduced) Aronsson system in fact consists 
of two linearly independent differential operators because of the perpendicularity 
between [Hp]+ and Hp; in fact, one may split Au = 0 to 


Hp(-,u, Du) D(H(-, u, Du)) 
H(-,u, Du) [Hp(-, u, Du)]- (Div (Hel, u, Du)) — H,(-, u, Du)) 


ll 
a) 


Theorem 5.1.4 makes clear that Aronsson’s absolute minimisers do not charac- 
terise the Aronsson system when N > 2, at least when the additional natural 
assumptions hold true. This owes to the fact that, unlike the scalar case, the 
Aronsson system admits arbitrarily smooth non-minimising solutions, even in the 
model case of the oo-Laplacian. For details we refer to [66]. 


Since Aronsson’s absolute minimisers do not characterise the Aronsson system, 
the natural question arises as to what is their PDE counterpart. The next theorem 
which is our second main result answers this question: 


5.1.5 Theorem [Divergence PDE characterisation of Abso- 
lute minimisers] 


Let wu: R" DQ — RN be a map in C1(Q; RY). Fix also O € 2 and consider the 
following statements: 


(I) u is a vectorial minimiser of E,.(-,O) in C!(O;RY)'. 


‘We remind the reader that u € Cj(O;R%) means u— g € Cj(O;R). 
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(II) We have 


max |Hp(-,u, Du) : Dy + H,(-,u, Du) - 2} > 0, 


Argmax{H(-,u,Du) : O} 


for any ~ € Ci(O;RY). 


(III) For any = € Cj(O; RY), there exists a non-empty compact set 


Ky, =K C Argmax{H(-,u, Du) : O} (5.1.10) 


such that, 


(Hp(-,u, Du) Db + H,(-,u, Du) wy) = 0. (5.1.11) 
Then, (I) => (I) => (III). If additionally H(z,-,-) is convex on RY x R‘*” 
for any fixed x € Q, then (III) => (I) and all three statements are equivalent. 
Further, any of the statements above are deducible from the statement: 


(IV) For any Radon probability measure o € P(Q) satisfying 
supp(a) C Argmax{H(-, u, Du) ‘ O}, (Oat. 12) 


we have 
— div(Hp(-,u, Du)o) + H,(-,u, Du)o = 0, (5.1.13) 


in the dual space (Ci(O; RY))*. 


Finally, all statement are equivalent if K = Argmax{H(-,u,Du) : O} in (III) (this 
happens for instance when the argmax is a singleton set). 


The result above provides an interesting characterisation of Aronsson’s concept 
of Absolute minimisers in terms of divergence PDE systems with measures as 
parameters. The exact distributional meaning of (5.1.13) is 


i (HC, Du :Dy + H,,(-, u, Du) . do _ 0 


for all» € C}(O;R%), where the “:” notation in the PDE symbolises the Euclidean 
(Frobenius) inner product in RY*". 


The idea of Theorem 5.1.5 is inspired by the paper [40] of Evans and Yu, wherein 
a particular case of the divergence system is derived (in the special scalar case 
N = 1 for the oo-Laplacian and only for 2 = O), as well as by new developments 
on higher order Calculus of variations in L° in |[61, 64, 70}. 


Note that, it does not suffice to consider only Q = O as in [40] in order to 
describe absolute minimisers. For a subdomain O C Q, it may well happen that 
the only measure o “charging” the points of O where the energy density H(-, u, Du) 
is maximised is the Dirac measure at a single point x € 0O. This is for instance 
the case for the standard “Aronsson solution” of the oo-Laplacian on R?, given by 
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u(x,y) = |a|*/° — |y|*/’, as well as for any other oo-Harmonic function which is 
nowhere Eikonal (i.e. |Du| is non-constant on all open subsets). 


We conclude this introduction by noting that the two vectorial variational con- 
cepts we are considering herein (Definitions 5.1.1-5.1.2) do not exhaust the plethora 
variational concepts in L°. In particular, in the paper [76] the concept of tight 
maps was introduced in the case of H(z,7, P) = ||P|| where || - || is the operator 
norm on R%*", Additionally, in the papers [14, 56] a concept of special affine vari- 
ations was considered which also characterises the Aronsson system, in fact in the 
generality of merely locally Lipschitz D-solutions. Finally, in the paper [12] new 
concepts of absolute minimisers for constrained minimisation problems have been 
proposed, whilst results relevant to variational principles in L°° and applications 
appear in [21, 22, 26, 45, 73, 74]. 


5.2 Proofs and a maximum-minimum principle 
for H(., u, Du) 


In this section we prove our main results Theorems 5.1.4-5.1.5. Before delving into 
that, we establish a result of independent interest, which generalises a correspond- 
ing result from [52]. 


5.2.1 Proposition [Maximum-Minimum Principles| 


Suppose Let u € C?(Q; RY) be a solution to (5.1.8), such that H satisfies 
(a) Hp(-,u, Du) has full rank on Q, 
(b) there exists c > 0 such that 

(€"Hp(e,n, P)) - (€"P) 2 e€"He(e.9, Pl’, 


for all € € R% and all (z,n, P) € 2 x RX x RN". 


Then, for any O € 2 we have: 


sup H(-,u, Du) = max H(-,u, Du), (5.2.1) 
Oo 00 
inf H(-,u, Du) = min H(-,u, Du). (pide) 


The proof is based on the usage of the following flow with parameters: 
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5.2.2 Lemma 


Let u € C?(Q;R%). Consider the parametric ODE system 


V(t) a €'Hp(-,u, Du) ’ t#0, 
Li (5.2.3) 
(0) = 2, 
for given c € and € € R%. Then, we have 
d T 
= (H(-,u,Du)|,) = €THp(-,u, Du) D(H(-,u,Du)) |. «, (5.2.4) 
d 2 
af uly) > cle He(-,u, Du) (6.2.5) 


5.2.3 Proof of Lemma 5.2.2 


The identity (5.2.4) follows by a direct computation and (5.2.3). For the inequality 
(5.2.5), we have 


“eT u(r(t)) = (Du(vit))) -40 


The lemma ensues. 


5.2.4 Proof of Proposition 5.2.1 


Fix O € Q. Without loss of generality, we may suppose O is connected. Consider 
first the case where rk(Hp(-, u, Du)) =n < N. Then, the matrix-valued map 
Hp(-,u, Du) is pointwise left invertible. Therefore, by (5.1.8), 


(Hp(-,u, Du)) Hp(-, u, Du) D(H(-,u, Du)) = 0 


which, by the connectivity of O, gives H(-,u,Du) = const on O. The latter 
equality readily implies the desired conclusion. Consider now the case where 
rk(Hp(-,u, Du)) = N <n. Fix x € O and a unit vector € € R” and consider 
the parametric ODE system (5.2.3) of Lemma 5.2.2. By the fullness of the rank 
of Hp(-,u, Du)), we have that 


|€"Hp(-,u, Du))| > c >0 on O. 
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We will now show that the trajectory y(t) reaches 0O in finite time. To this end, 
we estimate 


d 
|Dultz=(oydiamn(O) > [Dulioylrte) ~ 200) = |Z 
for some t € (0,t), by the mean value theorem. Hence, 


= |€"™Du(y(t)) - +(0)| t 
= |e Duta) - (7 HeC-u,Da)] oo) |t 


> cE" Hp(-,u, Du 


2 
t 


Ne 
> (eer) é. 


This proves the desired claim. Further, since u solves (5.1.8), by (5.2.4) of Lemma 
5.2.2 it follows that H(-,u,Du) is constant along the trajectory. Thus, if x € O 
is chosen as a point realising either the maximum or the minimum in QO, then by 
moving along the trajectory, we reach a point y € 0O such that H(-, u, Du)|, = 
H(-, u, Du)|,.- This establishes both the maximum and minimum principle. The 
proposition ensues. 


5.2.5 Remark [Danskin’s theorem] 


The central ingredient in the proofs of Theorems 5.1.4-5.1.5 is the next consequence 
of Danskin’s theorem: for any O € (2 and any u,¢ € C'(Q;R), we have the 
identities 


FhicgEm let 14,0) = max (He(-,u, Du) : Dé + H,(-,u, Du) - 4), 
dt t=0+ O(u) Fe ; r 
d . . . . 
dt ae Pg) Olu) (Het. u, Du) : Dé + H,(-,u, Du) - 0), 

where 


O(u) := Argmax{H(-,u,Du) : O}. 
Indeed, by [34, Theorem 1, page 643] and the chain rule we have 


d d 
az ot noo lu + td, O) = ay 


ae (max H(., u+td,Du+ {Dé) ) 


H(-,u + tg, Du + 1D) ) 


t=0t 


= Max < 
~ Ou) \ dt 


= max (Hp(-,u,Du) : Dé + H,(-,u,Du) -@). 
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This establishes the first identity of (5.2.6). The second one follows through the 
substitutions @ ~~ —¢, t~ -t. 


Now we may establish Theorem 5.1.4. 


5.2.6 Proof of Theorem 5.1.4 


(I) Suppose first that u is a rank-one absolute minimiser on 2. The aim is to 
show that (5.1.8) is satisfied on 2. This conclusion in fact follows by the results 
n [49], but below we provide a new shorter proof. To this end, fix x € ( and 
p € (0, dist(x, 0Q)) and let O := B,(x). We fix also € € R% and choose 


o(y) := &(ly—a|? —p’). 


Then, ¢ € Co (B,(x);span{é]). By Remark 5.2.5 and our minimality assumption, 
the definition of one-sided derivatives yields 


ia 
di 


E.(u+td,O0) <0< “| Exo(u + td, ). (5.2.7) 
dt |t=o+ 


t=0- 


Hence, by (5.2.7), (5.2.6) and continuity there exists a point x, with |x, — 2| < p 
which lies in the argmax set 


(B,(z))(u) = Argmax{H(-,u,Du) : B,(x)} 


such that 
= 0. (5.2.8) 


Xp 


(Ho(-,u,Du) : Dé + Hy(-,u,Du) - 4) 


Therefore, 
al (2He(-,u, Du)|,, (%p —2z) + H,(-, u, Du)|,, (Ip —z/— ’’)) =, (5.29) 
If x, lies in the interior of B,(x), then it is an interior maximum and therefore 
D(H(-,u, Du)) |, = 0. 


This means that (5.1.8) is satisfied at x,. If x, lies on the boundary of B,(x), then 
this means that 


Vy €B,(x), we have H(-, u, Du)|, < H(-,u, Du)|, 


Pp 
The above can be rewritten as 


By(x) © (ep) = {H(-,u,Du) < H(-,u,Du)|,, }, 


and note also that «, € 0B,(x) NOH(z,). Hence, the sublevel set H(x,) satisfied 
an interior sphere condition at x,. If D(H(-, u, Du)) | = 0 then (5.1.8) is again 


x 
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satisfied at x,. If on the other hand 


D(H(-,u,Du))|, A 0 


then OH(«,) is a C? manifold near x, and the gradient above is the normal vector 
at the point x,. Due to the interior sphere condition, this implies that this is also 
the normal vector to the sphere 0B,(x) at x,. Thus, there exists A 4 0 such that 


t)—« = AD(H(,u, Du)) = (5.2.10) 
By inserting (5.2.10) into (5.2.9) and noting that |x, — «| = p, we infer that 


al (Hp(,u, Du)D(H(-,u, Du)) )| 20: 


Lp 


By dividing by 2A and letting p — 0, we deduce that (5.1.8) is satisfied at the 
arbitrary xz € Q. 


Conversely, suppose that wu satisfies (5.1.8) on Q, together with the additional 
assumptions of the statement. Fix O € 2 and ¢ € Cj (O;span|{é]). Without loss 
of generality, we may suppose O is connected. Since ¢ = (€'@¢)€, for convenience 
we set g := €'¢ and then we may write @ = g€ with g € Ci(O). Then, the 
matrix-valued map Hp(-, Du) is pointwise left invertible. Therefore, by (5.1.8) 


(Hp(-, Du))"Hp(-, Du) D(H(-, Du)) = 0 on O, 
which, by the connectivity of O, gives 


H(-, Du) = const on O. 


Since g € C!(R") with g = 0 on 0O, there exists at least one interior critical point 
x € O such that Dg(z) = 0. By the previous, we have 


Es (tO) = He, Du(z)) 
= H(z, Du(z) + € @ Dg(z)) 
= H(z, Du(z) + De(z)) 
< sup H(z, Du(a) + D¢(z)) 
= E,.(u+¢,0). 


The conclusion ensues. 


(II) Suppose that wu has co-minimal area. Fix x € 2 and p € (0, dist(xz,0Q)). Fix 
ig © N(He(, U, Du)"|,.); 


noting also that by Definition 5.1.3 the above set is the reduced nullspace of 
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Hp(-,u,Du)' at x. This implies that there exists a C! extension € € C1(R"; RY) 
such that €(x) = € and (€)'Hp(-,u,Du) = 0 on the closed ball B-(x) for some 
e € (0,p). By differentiating the relation (€)'Hp(-,u,Du) = 0 and taking its 
trace, we obtain 


€-div(Hp(-,u, Du)) + Dé: Hp(-,u, Du) = 0, (5.2.11) 


on B.(x). Since u has oo-minimal area and € is an admissible normal variation, 
by using Remark 5.2.5 and arguing as in the beginning of part (I), it follows that 


= 0 (5.2.12) 


Le 


(é- H,,(-, u, Du) + DE: Hp(-,u,Du)) 


for some x, € (B.(x))(u), where 
(B-(x))(u) = Argmax{H(-,u,Du) : B.(x)}. 
By (5.2.11)-(5.2.12), we infer that 


= 0 


Le 


E(x.) - (div (Hp(-,u, Du)) = H,(-,u, Du)) 


and by letting ¢ — 0, we deduce that 


= 0, 


x 


&- (div (Hp(-,u, Du)) — Hy,(-,u,Du)) 


for any € € N(Hp(-, U, Du)"|,). Hence, wu satisfies (5.1.9) at the arbitrary x € 2. 


Conversely, suppose that. u solves (5.1.9) on Q. Fix O € 0 and ¢ € C!(R"; RY) 
such that ¢'Hp(-,u, Du) = 0 on O. Note further that by the continuity up to the 
boundary of all functions involved, the latter identity in fact holds on O. By the 
satisfaction of (5.1.9) and Definition 5.1.3, it follows that 


de (div (Hp(-, u, Du)) = H,(-,u,Du)) =), 
on O CQ. By differentiating ¢'Hp(-, u, Du) = 0, we obtain 
o- div(Hp(-, u, Du)) + D@: Hp(-,u,Du) = 0, 
on O. By the above two identities, we deduce 
@-H,(-,u, Du) + Dé: Hp(-,u, Du) = 0, 


on O. Since O(u) C O, Remark 5.2.5 yields that u is a critical point since the 
left and right derivative of E(u + t¢,O) at t = 0 coincide and vanish. Since by 
assumption H(z,-,-) is convex on RY x R**", it follows that E.(-,O) is convex 
on C!(O;R). Hence, the critical point u is in fact a minimum point for this class 
of variations. This establishes our claim. 
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(III) This is an immediate corollary of items (I) and (II). 


Now we conclude by establishing Theorem 5.1.5. 


5.2.7 Proof of Theorem 5.1.5 


Fix O € 0 and u,¢ € C1(Q;R%). We show that (I) => (II) => (II) and that 
(III) => (1) under the additional convexity assumption. By recalling Remark 
5.2.5, note that if 


Eo(u+td,O) > En(u,O), for allt € R, (0.2.18) 


then directly by (5.2.13) and the definition of one-sided derivatives, we have 


d d 
=| 5 <0<—| Ey(uttd,O). 5.2.14 
Slcq Bot + 46,0) 505 S| Bu(u+td,0) (5.2.14) 


This shows (I) ==> (II). If (II) holds, note that one also has that 


min _ |Hp(-,u, Du) : Dd + H,(-, u, Du) - a = 0, 
Argmax{H(-,u,Du) : O} 


for any ¢ € Cj(O;R). By (5.2.6) we see that (5.2.14) is satisfied and by continuity 
we obtain the existence of a non-empty compact set K = Kg C O(u) such that 


(Hp(-,u, Du) :D¢é + H,(-,u, Du) - 6)| ai, (5.2.15) 


Hence, (III) ensues. If now (5.2.15) holds true for some non-empty compact set 
K C O(u), then by (5.2.6) we have that (5.2.14) is true. If further H(z, -,-) 
is convex for all x € 2, then by Lemma 5.2.8 given right after the proof, t > 
E(u + td, O) is minimised at t = 0 and (5.2.13) holds true. 


(IV) => (III): Let ¢ € P(O) be any Radon probability measure satisfying (5.1.12). 
Then, by assumption 


[ (Het Du): Dé + H,(-,u, Du) - ¢)do =i 


for all 6 € Ci(O; RY). Fix any point  € O(u). By choosing the Dirac measure 


a € P(O) given by 


C= Oz 
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which evidently satisfies supp(7) = {Z} C O(u), we obtain 


for any Z € O(u). The conclusion ensues with K = O(u). 


(Hp(-,u, Du) :Dé + H,(-,u,Du) - 0) . 


oO 
= 0, 


(III) => (IV): If we have K = O(u) and 


(He(-,u, Du): Dd + H,(-, u, Du) - )| 


K 


7 [ (Heb Du) Dé + Hy(-,u,Du) 6) do 


= 0, 


then for any Radon probability measure o € P(Q) with supp(c) C K, we have 


[ (Hen Du) : Dé + Hy(-,u, Du) - 6) do _ 


for all ¢ € Ci(O;R). Hence, we have shown that 


—div(Hp(-, u, Du)o) + H,(-,u, Dujo = 0, 


in the dual space (Cj(O; RY))*. 


The next result which was utilised in the proof of Theorem 5.1.5 completes our 
arguments. 


5.2.8 Lemma 


Let f : 


R — 


5.2.9 Proof of Lemma 5.2.8 


R be a convex function. If the one-sided derivatives f’(0*) exist 
and f’(0-) <0 < f’(0T), then f(0) is the global minimum of f on R. 


By the convexity of f on R, for any fixed s € R there exists a sub-differential 
ps € R such that 


f@)-—f(s) > p.(t—s), for allt e€ 


For the choice t = 0 and s > 0, we have 


R. (5.2.16) 


and note also that since convex functions are locally Lipschitz, the set (ps)o<s<1 
is bounded. Thus, since f’(0*) exists and is non-negative, the above inequality 
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yields 
0 < 7'(0") < mint p, < co. 
s—0t 


Hence, by passing to the limit as s > 0* in the inequality (5.2.16) for t > 0 fixed, 
we obtain f(t) — f(0) > 0. The case of t < 0 follows by arguing similarly. 
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Chapter 6 


Conclusions and future work 


6.1 Conclusions 


We would like to mention that this thesis is a collection of published papers 
presented as chapters consist of original results. This work includes new results in 
the field of Calculus of Variations in L°. The new results are improved previous 
theorems by generalising and relaxing some of the conditions. Chapter 2 and 
Chapter 5 are joint papers with my supervisor Dr. N. Katzourakis. Chapter 3 is 
a joint paper with Dr. N. Katzourakis and Dr. H. Abugirda. While chapter 4 is 
single author paper. 


The main result of Chapter 2 is that we characterise local minimiser of the 
following functional 


Eoo(u, O) := ess supH(-,u,Du), we Wer(O,R®), OE. 
O 


for appropriate classes of affine variations of the energy as generalised solutions of 
associated PDE system which plays the role of Euler-Lagrange equation. Similar 
result was proven for H(z,7,P) = |P|? in [53]. That makes our result a gener- 
alisation of result in [53] since the Hamiltonian function H depends not only on 
gradient of the function but also on the function itself and the domain. 


Chapter 3 is the joint paper with Dr. N. Katzourakis and Dr. H. Abugirda. The 
author of this thesis gave an idea which partly impacted on the proof of the main 
result “Rigidity and flatness of maps with tangential Laplacian in separated form” , 
which states let Q C R? be an open set and let also N > 2. Let u: R” DQ — RY 
be a classical solution to the nonlinear system 


[Du] Au=0 in Q, 


having the separated form u(z,y) = f(x) — f(y), for some curve f € (W??M 
C?)(R,R) and some p > 1. Then, the image u(Q) of the solution is contained in 
an at most countable union of affine planes in RY. 


67 


Chapter 4 investigated with authors own initiative and he found a new classi- 
cal oo-harmonic functions in high dimensions, particularly when domains are two 
dimensional in polar coordinates and at least two dimensional in Cartesian coor- 
dinates. The challenges were the technical computations and the regularity of the 
solutions on uncertain domains which was assumed to be well defined. 


We have two main outcomes of Chapter 5 which is a joint paper with Dr. N. 
Katzourakis. First result “Variational Structure of Aronsson’s system” coincides 
with the result in [49] when the Hamiltonian function depends only on the gradient 
function, namely H(x,7, P) = |P|?. The result characterises C? oo—minimal maps 
as solution of the Aronsson system and vice versa. One of the difficulty of proving 
this theorem was that we can not differentiate 


(€)'Hp(-,u, Du) = 0. 


if € ¢ C!(R"; RY). We avoided that using reduced nullspace in the definition 5.1.3. 
The second result “Divergence PDE characterisation of Absolute minimisers” is a 
completely new original result. Lets highlight the main differences of this theorem 
with previous results. Firstly it has been proved for C! maps and for C! variations 
vanishing on compactly contained boundaries. Secondly we did not use approxi- 
mation techniques of the L? space but rather techniques of the L° space, namely 
Danskin’s theorem. 


6.2 Future work 


We believe that the work in this field is interesting and there are still many open 
problems one can work on, for example: 


1. It is common that a solutions to a PDEs might have less regularity than we 
require. So it is natural to work on extending the result of the theorem 5.1.5 
from C! to Sobolev spaces or D-solutions. 


2. Theorem 5.1.5 gives us characterization of Absolute minimiser only on Argmax 
set. One of the methods to fill the gap is to study vectorial L°°-absolute 
minimisers on ( using vectorial L?-absolute minimisers on Q, i.e. for ev- 
ery O € 2°. if uy is minimiser of E,(u,O) := |H(-,u,Du)|z(0), then one 
can study “convergence” of uy tO Us as p — Oo, where wu. is minimiser of 
Exo (u, O) = |H(., U, Du) |r): 


3. All known explicit solutions have at least C'® regularity. So it is challenging 
to find explicit solutions of the theorem 5.1.5 and find out how and/or why 
other results have to fail. 
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